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CHAPTER  ONE 


INTRODUCTION 


1.1  Physical  Problem 

We  are  concerned  with  the  problem  of  passive  sonar  processing  in 
mobile  sonars,  sucn  as  towed  line  arrays.  The  physical  mechanism 
motivating  the  study  is  hydrodynamic  noise,  which  often  interferes 
with  sonar  (Sound  Navigation  Ranging)  performance.  It  defines  a  class 
of  noises  arising  from  the  flow  of  water  over  hydrophones  or  the 
hydrophone  housing  [1].  Flow  noise,  considered  a  member  of  this 
class,  remains  after  all  other  sources  of  hydrodynamic  noise  have  been 
accounted  for.  Its  energy  levels  increase  with  tow  speed,  and 
hydrophones  are  usually  directly  coupled  to  the  excitation.  Hence, 
flow  noise  is  often  a  major  source  of  interference  at  low  frequencies 
and  high  tow  vessel  speeds. 

The  flow  noise  components  of  primary  concern  to  us  are  caused  by 
turbulent  flow  in  a  turbulent  boundary  layer  (TBL).  The  TBL  existing 
between  tne  free  stream  flow  and  tne  array  creates  random  pressure 
fluctuations  and  transmits  tnem  to  the  hydrophones.  These  pressure 
fluctuations  are  usually  seen  as  additive,  often  predominant,  noise 
components,  independent  of  the  signal. 
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In  the  usual  sonar  application,  the  hydrophones  in  Fig.  l.l 
simultaneously  sense  the  signal  and  noise  pressure  field.  Hydrophone 
data  are  then  sent  to  data  processing.  The  processed  data  are 
presented  to  an  ooserver,  who  must  decide  whether  or  not  a  target  is 
present.  Uncertainty  in  this  decision  process  is  obviously  increased 
by  the  presence  of  corrupting  flow  noise.  The  problem  is  to  eliminate 
as  much  of  the  unwanted  flow  noise  as  possible  without  destroying 
target  information. 

Usual"  methods  of  reducing  direct-coupled  TBL  noise  (flow  noise) 
in  towed  arrays  fall  into  the  category  of  "mechanical"  processing. 

One  such  method  tnat  has  had  some  success  is  referred  to  as  hosewall 
attenuation.  It  embeds  tne  sensors  in  material  tnat  tends  to  flatten 
the  flow  noise  wavenumber  spectrum.  If  the  nosewall  sufficiently 
whitens  the  noise  spectrum,  then  conventional  oeamforming  on  tne 
sensor  outputs  is  optimum.  The  hosewall  may  also  sufficiently 
attenuate  the  iiign  wavenumber  components  so  tnat  flow  noise  is  no 
longer  the  predominant  noise  component.  When  the  remaining  noise  is 
spatially  white,  then  conventional  oeamforming  is  again  nearly 
optimal.  If  it  is  desired  to  improve  signal-to-noise  ratio,  then 
hosewall  attenuation  must  preferential ly  attenuate  flow  noise  while 
passing  the  signal  undistorted. 

We  shall  argue  on  the  oasis  of  an  assumed,  but  physically 
reasonable,  flow  noise  model  that  hosewall  attenuation  gains  nothing 
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over  conventional  processing  on  space-continuous  apertures  if  the 
aperture  contains  a  sufficient  number  of  noise  correlation  lengths. 
Thus,  space-continuous  extended  sensors  essentially  whiten  the  flow 
wavenumber  spectrum  with  an  insignificant  effect  on  signal.  Hence,  a 
detector  that  uses  conventional  beamforming  on  extended  sensor  outputs 
can,  with  noise  uncorrelated  from  sensor  to  sensor,  come  close  to 
achieving  the  best  possible  detection  index  with  an  extremely  simple 
but  "robust"  instrumentation.  The  present  investigation  grew  out  of 
earlier  work  by  a  numoer  of  researchers.  Pasupathy  [2]  analyzed 
space-time  processing  on  a  space-continuous  finite  observation 
interval.  He  established  asymptotic  conditions  on  noise  spatial 
parameters  under  which  optimal  processing  could  be  expected  to  yield 
significant  detection  improvement  over  simpler  conventional 
processing,  such  as  conventional  oeamforming.  In  the  same  manner,  he 
established  asymptotic  conditions  for  which  conventional  processing  is 
i tse If  nearly  optimal . 

Although  Pasupathy' s  results  serve  mainly  as  asymptotic  check 
points  for  parameter  values,  such  as  the  ratio  of  interference  power 
to  the  white  noise  power  in  the  interference  band  and  the  noise 
spatial  bandwidth  array  length  product  (sL),  they  do  not  provide 
sufficient  insight  into  the  practical  issues  of  tne  array  design 
problem.  Pasupathy's  model  included  an  additive  spatially  white  noise 
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component,  incoherent  from  point  to  point.  This  formally  results  in 
infinite  power  at  a  point;  hence,  it  is  not  a  realistic  model  of 
underlying  flow  noise  processes.  Nor  can  we  derive  meaningful 
inferences  from  Pasupatny's  results  applicable  to  the  detection 
problem  with  discrete  observations  (arrays). 

In  recent  years  there  nave  oeen  extensive  investigations  of  the 
TBL  and  its  effect  on  sonar  processing.  Much  of  the  work  is 
classified.  The  following  summarizes  the  most  relevant  unclassified 
work.  Jorgensen  and  Maidanick  [3]  have  examined  signal-to-noise 
improvement  in  tne  TBL  that  is  achievable  by  adjusting  sensitivities 
and  separations  of  two-  and  tnree-point  transducer  systems.  Kennedy 
[4]  considered  the  approach  of  Nuttall  and  Hyde  [5]  of  fully  optimal 
weighting  of  discrete  sensor  outputs  clustered  in  TBL  and  isotropic 
noise.  Motivated  by  these  earlier  studies.  Fay  and  Owsley  [6] 
proposed  an  extremely  simple  instrumentation  that  combined  the  outputs 
of  two  sensors,  resulting  in  a  simple  noise  cancellation 
instrumentation.  They  considered  the  concept  of  minimum  variance 
complex  weighting  on  clusters  of  point  sensors  for  the  class  of  TBL 
noise  fields  represented  oy  a  very  simple  exponential  model  with 
additive  uncorrelated  noise. 
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However,  a  number  of  issues  remain  unresolved.  For  instance, 
there  is  not  yet  an  adequate  characterization  of  the  achievable  gains 
to  be  made  with  optimal  processing  over  conventional  processing  in  a 
flow  noise  environment,  especially  with  regard  to  wideband  detector 
performance.  It  is  also  not  clear  wnat  detector  instrumentations  are 
necessary  to  achieve  the  gains.  Attainable  optimal  and  conventional 
performance  values  are  not  known  nor  are  tne  exact  conditions  on 
spectral  and  array  parameters  for  wnicn  optimal  processing  can  produce 
large  gains  over  conventional  processing.  Insight  into  these  issues 
is  necessary  for  practical  interpretation  of  proposed  signal 
processing  schemes  for  flow-noise-dominated  environments. 

Our  results  differ  from  previous  work  in  several  ways.  In 
formulating  the  detection  problem,  we  derive  exact  closed-form 
expressions  for  the  optimal  detector  structure  and  also  the  best 
detection  performance  achievaole  on  a  space-continuous,  finite 
observation  interval.  This  is  done  by  modeling  TBL  noise  with  a 
first-order  Butterworth  rational  wavenumber  spectrum  consistent  with 
TBL  models  used  in  [3],  [4],  and  loJ  and  also  with  Pasupathy's  model. 
Altnough  we  confine  our  attention  to  the  TBL  flow  noise  problem  in 
sonar  line  arrays  (towed  line  arrays),  the  general  results  have  a 
wider  application  to  other  interfering  noise  sources.  We  omit 
discussion  of  such  noise  sources  oecause  their  influence  on  sonar 
performance  nas  been  treated  extensively  in  the  literature  [7]  - 
[18].  We  shall  present  a  unified  treatment  of  optimal  and 
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conventional  detector  performance  in  passage  from  a  discrete  array  of 
spatial  measurements  to  a  space-continuous  observation  interval.  It 
is  shown  that  an  uncorrelated  noise  component  has  zero  spectral  level 
in  the  limit  of  continuous  observations.  Therefore,  only  spatially 
colored  flow  noise  remains  to  limit  detection  capability  on  the  space- 
continuous  interval.  Since  even  tne  slightest  amount  of  uncorrelated 
noise  brings  conventional  processing  closer  to  optimal,  our 
performance  results  represent  tfie  best  possible  performance  and  also 
the  absolute  upper-bound  performance  advantage  attainable  with  optimal 
processing  over  conventional  processing. 

Once  specific  spectra  have  been  assigned  to  tne  signal  and  noise 
processes,  the  binary  detection  problem  of  a  coherent  wavefront 
signal  in  spatially  colored  flow  noise  with  continuous  observations 
becomes  a  straightforward  filtering  problem.  The  object  then  is  to 
secure  maximum  signal-to-noise  ratio  subject  to  the  spatial  bandwidth 
constraint  imposed  Dy  tne  finite  length  array  at  each  signal 
frequency.  Determining  the  optimal  filter  (detector)  involves  solving 
Fredholm  integral  equations  encountered  in  connection  with  the 
likelihood  ratio  test  [19].  Tne  resulting  optimal  instrumentation  is 
realized  as  an  inner  product  between  the  space-continuous  data  and  the 
solution  of  tne  integral  equation. 
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This  optimum  instrumentation  obviously  depends  on  the  particular 
features  of  the  typical  flow  noise  spectrum  and  the  way  these  features 
interact  with  array  properties.  It  is  well  known  that  a  narrowband 
conventional  detector  (i.e.,  a  conventional  beamformer)  is  formally 
optimal  when  the  noise  wavenumber  spectrum  is  white.  We  shall  show 
that  array  gain  for  the  conventional  beamformer  is  nearly  optimal  when 
the  array  length  contains  a  large  number  of  noise  correlation  lengths. 
Thus,  the  long  array  must  reduce  the  problem  to  near-white  noise  in 
the  wavenumber  interval  for  which  the  beamformer  response  is 
significant. 


The  main  interest  here  is,  of  course,  in  situations  where  the 
noise  is  not  white  but  instead  has  a  general  bandpass  behavior  in 
wavenumoer  space.  In  such  situations,  gains  witn  optimal  processing 
over  conventional  processing  on  a  finite  interval  can  be  significant. 

Our  formulation  of  the  detection  problem  of  a  signal  in  flow 
noise  with  continuous  observations  reduces  to  a  classical  binary 
detection  problem  of  a  known  spatial  signal  in  Gaussi an-coloreo 
noise.  Location  of  the  signal  is  assumed  to  be  known.  Distance  to 
the  signal  from  a  reference  point  in  tne  array  is  assumed  to  be 
sufficiently  large  so  that  the  signal  is  in  tne  farfield  and  thus 
presents  a  plane  wave  pressure  field  to  the  array.  The  signal  is 
assumed  to  be  Gaussian  zero  mean  with  known  spectral  properties. 

Noise  is  assumed  to  be  Gaussian  zero  mean  with  known  temporal  and 
spatial  spectra.  Numerous  detection  problems  of  this  kind  have  been 
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solved  in  the  framework  of  temporal  processing  on  a  finite  temporal 
observation  interval  [19],  [20].  Here  we  deal  with  a  space-time 
interval  but  use  established  mathematical  procedures  to  solve  the 
detection  problem. 

The  binary  detection  proolem  wiui  arrays  of  equally  spaced 
discrete  sensors  encounters  an  unavoidable  aliasing  proolem  in  a  flow 
noise  environment.  Since  the  upper  cutoff  wavenumber  of  the  typical 
flow  noise  wavenumber  spectrum  is  very  high,  there  will  be  severe 
folding  of  noise  wavenumoers  onto  acoustic  wavenumbers  unless  uniform 
sensor  spacing  is  extremely  tight.  The  resulting  aliased  noise  level, 
at  tne  signal  wavenumber,  is  thus  increased  over  the  noise  level  with 
continuous  observations;  also,  the  aliased  noise  wavenumber  spectrum 
is  often  flatter.  As  a  consequence,  optimal  and  conventional 
detection  capabilities  are  degraded  substantially  relative  to  the  best 
continuous  performance.  With  the  flatter  spectrum,  tne  conventional 
processor  itself  moves  closer  to  being  optimum.  The  spatial  Nyquist 
rate  of  the  inverse  noise  spatial  bandwidth  [21]  associated  with 
sampling  a  bandpass  flow  noise  process  is  meaningless  in  the  spatial 
case  since  it  implies  the  generation  of  a  Hilbert  transform  on  a 
space-continuous  interval.  Once  discrete  sensors  are  introduced,  the 
aliasing  loss  has  happened  and  is  irreversible. 

In  almost  all  sonar  applications,  the  number  of  available  sensors 
or  telemetry  channels  is  limited  and  not  sufficient  to  fully  populate 
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a  given  observation  interval.  This  raises  the  question  of  efficient 
use  of  a  given  number  of  sensors  aimed  at  achieving  the  best  possible 
detection  performance.  We  show  that  there  are  better  ways  of  using  a 
limited  number  of  sensors  than  simply  spacing  them  uniformly,  making 
the  noise  uncorrelated  from  sensor  to  sensor.  We  consider  the  array 
construction  method  in  [6],  consisting  of  groups  of  tightly  spaced 
sensors  with  noise  correlated  between  sensors.  By  using  such  groups, 
one  can  get  a  much  higher  array  gain  at  certain  frequencies  and  hence 
achieve  greatly  improved  detector  performance  (either  optimal  or 
conventional)  over  the  performance  acnievable  with  the  same  number  of 
sensors  spread  out  in  the  observation  interval  so  as  to  make  the  noise 
uncorrelated  from  sensor  to  sensor. 

Finally,  we  note  that  array  gain  is  frequently  used  to  compare 
the  performance  of  alternative  array  processing  procedures.  It 
specifies  the  postbeamforming  signal-to-noise  ratio  at  any  given 
frequency  in  terms  of  the  input  signal-to-noise  ratio  and  is, 
tnerefore,  independent  of  the  transfer  function  of  a  linear  filter 
connected  to  the  beamformer  output.  Wideband  detector  or  tracker 
performance,  on  the  other  hand,  is  not  independent  of  such  a  filter. 
Thus,  even  if  two  systems  have  the  same  array  gain,  their  wideband 
performance  can  De  quite  different.  Specifically,  in  the  case  of 
towed  line  arrays  operating  in  flow  noise  tnat  is  heavily  concentrated 
at  very  high  wavenumbers,  optimal  and  conventional  detectors  can  have 
very  different  wideband  detection  indices.  It  is  shown  here  that  much 
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of  tne  performance  loss  due  to  the  use  of  conventional  beamforming  can 
be  recovered  by  a  simple  modification  of  the  Eckart  filter  usually 
placed  at  the  beamformer  output.  For  a  uniform  array  that  is  long 
compared  with  the  flow  noise  correlation  distance,  such  a  filter  can 
reduce  the  difference  between  optimal  and  conventional  beamformer 
performance  to  a  negligible  amount.  Tne  modified  Eckart  filter  should 
be  much  easier  to  implement  tnan  a  fully  optimal  Deamformer  and,  at 
least  in  environments  such  as  those  considered  here,  should  realize 
much  of  the  available  performance  gain. 

It  is  appropriate  at  this  point  to  discuss  briefly  the  general 
form  of  the  flow  noise  model  to  be  used  in  this  study,  saving  the 
detailed  analytic  description  for  Chapter  2.  It  is  known  that  the  TBL 
pressure  fluctuations  are  random  in  both  space  and  time,  but  the 
underlying  physical  mechanisms  are  not  well  understood  [22,  p.  125]. 
This  lack  of  complete  knowledge  is  not  overly  restrictive  in  our 
case.  We  need  only  consider  available  measured  statistical  properties 
of  the  TBL  processes  in  order  to  specify  the  optimal  processor. 

We  shall  characterize  the  TBL  process  Dy  a  zero  mean,  Gaussian 
model,  completely  described  by  its  covariance  function.  There  is 
reason  to  believe  that  TBL  pressure  fluctuations  are  Gaussian  [22, 
p.  125].  More  pragmatically,  tne  Gaussian  property  greatly  simplifies 
analysis,  leading  to  a  linear  processor  that  is  optimum  for  a  large 
class  of  error  criteria  [19]. 
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One  may  obtain  insight  into  tne  TBL  pressure  fluctuations  from  a 
discussion  by  Strasberg  [22].  Hydrophones  (pressure  sensors)  in  the 
moving  array  (see  Fig.  1.1)  respond  to  the  local  random  pressure 
fluctuations  generated  by  the  turbulence.  In  the  case  of  towed 
arrays,  the  turbulent  flow  results  mainly  from  the  motion  of  the  tow 
vessel  and  hence  from  the  array  moving  tnrough  the  water.  Turbulent 
motion  is  characterized  by  flow  velocities,  whose  magnitude  and 
direction  vary  randomly  witn  time.  These  flow  velocities  are 
superimposed  on  the  mean  flow  velocity,  which,  for  the  towed  array,  is 
of  the  order  of  a  typical  tow  vessel  speed  (17-42  ft/sec). 

Local  velocity  variations  of  the  eddies  associated  with  the 
turbulence,  traveling  at  the  mean  flow  velocity  along  the  array,  are 
sensed  at  various  temporal  frequencies.  If  Uc  is  the  mean 
(convective)  speed  of  the  moving  sensor  relative  to  the  flow,  the 
wavenumber,  v.  ,  associated  with  frequency  u>  =  2 -nf.  is 

K  K 


These  local  velocity  variations  give  rise  to  pressure  fluctuations 
along  the  array  at  a  frequency  f^.  Tnese  pressure  fluctuations  are 
associated  with  alternating  hign  and  low  pressure  values  along  the 
array.  The  correlation  function,  describing  the  velocity  and  hence 
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the  pressure  dependence  between  points  in  the  array  at  frequency  f^, 
has  an  oscillatory  behavior  witn  spatial  period  x^  =  Since 

typical  mean  flow  speed  is  of  tne  order  of  25  ft/sec,  adjacent  zeros 
of  the  correlation  function  are  separated  by  distances  of  the  order  of 
i nches . 

The  local  velocity  variations  of  tne  eddies  are  commonly  believed 
to  maintain  their  spatial  shape  or  structure  for  some  distance  as  they 
move  along  the  array  at  the  mean  flow  speed  in  the  direction  of  flow. 
However,  they  eventually  lose  their  structure  and  thus  become 
uncorrelated  over  distances  that  are  often  very  large  compared  with 
the  correlation  period  \  but  small  compared  with  acoustical  signal 
wavelengths  at  typical  sonar  frequencies. 

Theoretical  models  of  TBL  spectra  do  exist  [23],  [24,  p.  145]  but 
do  not  agree  well  with  measured  spectra  for  general  towed  array 
conditions.  Experimental  measurements  [25],  [26j  suggest  that  tne  TBL 
wavenumber  spectrum  assumes  tne  form  shown  in  Fig.  1.2.  This 
representation  also  agrees  witn  the  experimental  Corcos  spectral  model 
for  TBL  wall  pressure  [27],  [28].  Available  measurements  assure  us 
that  the  typical  wavenumber  spectrum  is  concentrated  about  the 
convective  wavenumber 
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unce  again,  convective  speed  is  of  the  order  of  the  free  stream  speed 

(e.g.,  tow  vehicle  speed),  which,  for  sonar  applications,  is  orders  of 

magnitude  slower  than  the  speed  of  sound  in  water.  Thus,  v  is 

k 

likely  to  be  a  high  wavenumber,  far  removed  from  the  acoustic  region 
where  signal  and  other  noise  components  with  wavenumbers  characterized 
by  the  speed  of  sound  are  found.  It  is  this  large  wavenumber 
difference  between  signal  and  flow  noise  that  provides  the  principal 
mechanism  for  separating  signal  and  interfering  flow  noise  and,  hence, 
potentially  achieving  large  detection  gains.  The  peak  noise  level  at 
vk  is  typically  orders  of  magnitude  larger  than  the  noise  level  at 
the  signal  wavenumber.  The  wavenumber  spread  about  the  peak  is 
character i zed  here  by  the  spatial  bandwidth,  sk,  which  characteri zes 
the  correlation  scale  of  the  TBL  pressure  fluctuations  along  the 
array.  Experimental  results  in  [26]  suggest  tnat  6k  nas  the  form 


Since  both  vk  and  bk  vary  linearly  with  ,i>  ,  the  fractional 
Dandwidth  a  is  frequency  independent.  The  convective  wavelength  is 
V  which  is  frequency  dependent.  Experimental  results  [26] 
indicate  that  components  of  pressure  fluctuations  in  the  turbulent 
flow  at  frequency  wk  are  essentially  uncorrelated  after  traveling 
five  convective  wavelengths  in  the  direction  of  flow.  This 
corresponds  to  values  of  a  of  order  0.02. 
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The  small  correlation  lengtn  scale  of  the  turbulent  pressure 
fluctuations,  shown  in  Fig.  1.2c,  has  important  practical  implications 
in  sonar  applications  where  one  tends  to  separate  sensors  by  one-half 
of  an  acoustical  wavelength.  In  such  cases,  one  sees  that  the  high 
wavenumber  TBL  noise  process  is  essentially  uncorrelated  from  sensor 
to  sensor.  This  generally  results  in  aliased  flow  noise  wavenumber 
density  everywhere  in  the  acoustic  wavenumber  interval.  Hence, 
detection  capability  is  severely  degraded  relative  to  detection 
capaDility  with  extremely  tight  sensor  spacing. 

It  is  clear,  as  we  shall  demonstrate,  that  there  are  preferred 
sensor  spacings  for  the  typical  flow  noise  spectrum  in  Fig.  1.2.  For 
instance,  if  the  sensor  interval  in  Fig.  1.2c  is  Wv^,  noise  tends 
to  cancel.  Sensor  spacings  of  2W are  bad  because  the  flow  noise 
adds  coherently.  This  suggests  that  sensors  should  always  be  tigntly 
spaced  within  the  noise  correlation  length,  Dut  that  the  precise 
spacing  is  not  unimportant  if  one  desires  to  improve  detection 
perf ormance  in  a  flow  noise  environment. 

The  flow  noise  power  spectrum  is  thougnt  to  be  flat  at  low 

-3 

temporal  frequencies  but  to  slope  at  a  rate  of  approximately  f  at 
high  frequencies.  The  transition  between  low  and  high  frequencies, 
given  by  the  quotient  of  the  free  stream  velocity  to  the  Doundary 
layer  thickness,  is  usually  in  the  hundreds  of  hertz  [1,  pp. 

360-365].  We  shall  assume  that  the  flat  portion  of  the  power  spectrum 
is  either  coincident  with  or  larger  than  the  signal  band  of  interest. 
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Thus,  we  are  always  dealing  with  a  flat  noise  spectrum  in  the  signal 
band  with  a  power  level  that  varies,  according  to  [1],  as  the  cube  of 
the  free  stream  speed. 

1.2  Organization 

The  study  covers  material  tnat  can  be  divided  into  six  parts. 
Chapter  2  provides  the  mathematical  background,  models,  and 
terminology  needed  in  the  remaining  chapters. 

Chapter  3  examines  the  binary  detection  problem  of  a  random 
conerent  wavefront  signal  in  Gaussian  stationary  flow  noise  modeled  as 
a  first-order  Butterworth  rational  wavenumber  spectrum.  We  evaluate, 
in  detail,  the  performance  and  structure  of  the  lixelihood  ratio  test 
by  solving  Fredholm  integral  equations  over  a  finite  and  continuous 
observation  interval.  This  characterizes  the  spatial  aspects  of  the 
optimum  detector  and  provides  definite  values  for  the  detection  index 
in  terms  of  tne  spectral  parameters  of  signal  and  flow  noise  and  their 
interaction  with  the  aperture  parameters. 

In  Chapter  4  the  continuous  observation  interval  is  discretized 
to  model  an  array  of  equally  spaced  point  hydrophones.  Again,  we 
specify  the  optimal  instrumentation .  Array  gain  and  detection  index 
are  used  to  measure  comparative  detection  performance  of  optimal  and 
conventional  detectors  in  tne  flow  noise  environment.  It  is  shown 
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that  much  of  the  performance  loss  due  to  the  use  of  conventional 
beamforming  can  be  recovered  by  a  simple  modification  of  the  Eckart 
filter  usually  placed  at  the  beainformer  output. 

Chapter  5  deals  with  array  design  considerations,  notably  the 
choice  of  array  geometry.  When  sensors  are  continuous,  the  primary 
constraints  concern  their  length  and  separation.  When  sensors  are 
discrete,  one  is  almost  certainly  restricted  to  a  given  number  of 
sensors  and  must  arrange  these  in  a  geometrical  pattern,  which 
maximizes  their  detection  index,  preferably  without  unduly 
complicating  the  instrumentation.  Both  optimal  solutions  and  the 
compromises  dictated  by  practical  considerations  are  presented. 

Chapter  6  summarizes  the  results  and  recommends  research  along 
related  lines. 
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ORIGIN 


Fig.  1.1.  Linear  Array  With  Plane  Wave  Signal  and  Noises 
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CHAPTER  TWO 


MATHEMATICAL  BACKGROUND 


2.1  Mathematical  Model 

This  cnapter  develops  quantitative  signal  and  noise  models  with 
mathematical  properties  to  be  used  throughout  the  dissertation.  It 
formulates  the  detection  problem  for  discrete  and  continuous  spatial 
observations  and  discusses  the  resulting  integral  equations. 

The  space-time  processing  model  for  target  (signal)  in  a  flow 
noise  environment  is  shown  in  Fig.  2.1.  The  received  pressure  field 
is  observed  over  the  entire  space-continuous  array.  Each  point,  x,  in 
the  array  receives  a  space-time  waveform, 

y(t,x)  =  n(t,x)  +  A  s ( t,x ) ,  (2.1) 

L  L 

~  2  -  x  -  2 

T  *  T 

2  <  t  ^  5 

wnich  consists  of  signal  s(t,x)  and  additive  noise  n(t,x).  A  is  1  if 
signal  is  present  and  0  if  it  is  absent.  The  received  waveform  is  a 


realization  of  a  random  process  parameter i zed  by  real  numbers  t  and 
x.  The  realization  contains  target  information  when  a  target  is 
present,  and  the  noise  usually  consists  of  several  components,  such  as 
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uncorrelated  noise  (electronic  noise),  farfield  i nterference,  and  flow 
noise  induced  by  movement  of  water  past  a  mobile  sonar  array.  We  take 
special  interest  in  the  flow  noise  problem  associated  with  high-speed 
towed  line  arrays  and  will  generally  ignore  other  noise  components 
whose  influence  on  the  detection  proDlem  has  been  studied  extensively 
[7]  -  [18]. 

2.2  Data  Representation 

2.2.1  Fourier  Series  Expansion 

The  received  waveform  is  ODserved  over  a  temporal -spati al 
observation  interval,  [-T/2,  T/2]  x  [-L/2,  L/2].  The  usual  question 
in  detection  is  whether  or  not  the  target  signal  is  present  in 
y(t,x).  This  is  a  binary  decision  proolem,  with  the  following 
hypotheses  on  the  received  waveform  (Van  Trees  [19]): 


H,  :  y(t,x)  =  n(t,x)  +  s(t,x)  ^Target  present) 

H  :  y(t,x)  =  n(t,x)  (Noise  only  present)  . 

0 

Signal  and  noise  in  (2.2)  will  oe  modeled  as  zero  mean  Gaussian 
random  processes  statistically  independent  of  eacn  other . 
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A  data  vector,  jr,  is  constructed  by  representing  the  received 

space-time  pressure  field  in  any  convenient  manner.  An  obvious, 
though  not  necessarily  desiraole,  choice  is  a  succession  of  space-time 
samples,  yU^x.): 


yUpXj) 

y ( t i ’ xm ) 

y( tN’X1) 

y(tN5xM) 


(2.3) 


where  i  =1,  2,  ...,  N  and  j  =  1,  2,  . . . ,  M  so  that  there  are  N  time 

samples  at  each  of  M  spatial  measurement  points.  The  components  of  r 

are  zero  mean  and  jointly  Gaussian  so  tnat  their  statistical 
properties  are  specified  by  a  covariance  matrix  with  components: 


y ( t . ,x, ) 


★ 

y  U, 


’ xm^ 


(2.4) 


Tne  representation  in  u\3)  and  (2.4)  is  cumbersome  since  it 
requires  closely  spaced  time  samples.  For  an  observation  interval  T 
and  bandwidtn  w,  one  generates  a  covariance  matrix  of  dimension  (M  x 
2wT)  witn  little  useful  structure.  The  syinool  *  denotes  a  transpose 
conjugate  operation  and  E  represents  statistical  expectation. 
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A  more  convenient  representation  is  available  when  all  of  the 
following  conditions  are  satisfied:  (a)  the  Tw  product  of  the  signal 
is  large,  a  condition  usually  associated  witn  passive  sonar  detection 
problems,  (b)  the  signal  wavefront  travels  across  tne  aperture  (array) 
in  a  time  that  is  snort  compared  with  T,  and  (c)  the  turbulent 
boundary  layer  noise  travels  its  correlation  distance  in  a  time  that 
is  short  compared  with  T  and  tne  noise  correlation  time  is  also  short 
compared  with  T.  Tnen,  Fourier  coefficients  at  different  frequencies 
are  approximately  uncorrelatea  [29j,  [30].  Tne  Fourier 
representation ,  therefore,  generates  a  block  diagonal  covariance 
matrix  in  (2.4)  and  nence  reduces  the  effective  dimensionality  to  M. 

The  data  vector  is  now 


r 


(2.5) 


where 


dt 


2irk 


(2.6) 


k  =  1,  2,  ...  N  . 


22 


TR  7997 


The  received  waveforms  are  real  so  tnat  Y(-u>.  ,x.)  =  Y*(<d.  ,x.)«  If 

K  J  K  J 


N  _>  2Tw,  then  r.  should  be  an  adequate  representation  of  the  continuous 
data  over  (-T/2,  T/2).  The  elements  of  r  are  jointly  Gaussian  since 
the  Fourier  series  representation  corresponds  to  a  linear 
transformation  on  time  samples.  The  process  is  thus  characterized  by 
its  (NM)x(NM)  covariance  matrix: 


A 


E  [r  r  ]  . 


K 


(2.7) 


Since  Fourier  coefficients  associated  with  different  frequencies 
are  uncorrel ated. 


K 


(2.8) 


n,m 


(2.9) 


(2.10) 
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and 


K(u  )  = 

v  rr 


^I(«n)iun)*j  ,  n  = 


1,  2, 


(2.11) 


Consequently,  the  covariance  matrix  in  (2.7)  nas  the  following  block 
diagonal  form: 


o  .  .  .0 


o  .  .  .  o  K. 


(2.12) 


where  block  matrices  are  defined  by  K .  =  K(w.)  =  K(2nj/T). 

J  3 


2.3  Target  Model 

The  target  (signal)  component  of  tne  received  waveform  is  modeled 
as  shown  in  Fig.  2.2.  Ttie  signal,  s(t),  originates  at  a  source  point 
that  is  assumed  to  satisfy  farfield  conditions.  Tne  signal  pressure 
field  propagates  as  a  coherent  wavefront  from  the  source  point,  p,  in 
the  direction  of  the  unit  vector,  5_,  making  an  angle,  9  ,  with  the  x 
axis,  which  coincides  with  the  coordinate  axis  for  the  array.  The 
following  properties  are  assumed  for  the  signal: 
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1.  The  signal  s(t)  is  a  stationary  zero  mean  Gaussian  random 
process  characteri zed  oy  known  second  moment  properties. 

2.  Tne  signal  s(t)  propagates  undistorted  from  p  to  x,  where  it 
is  represented  as  s(t,x),  x  e(-oo,a >). 

3.  The  signal  s(t,x)  is  statistically  independent  of  the  noise 
processes  n(t,x). 

4.  T  is  large  compared  with  tne  signal  correlation  time  and 
signal  travel  time  across  tne  spatial  ooservation  interval 
(array). 

Fourier  coefficients  of  the  signal  are  given  by 


S(u>k,x) 


A  •  A 

where  s(t,x)  =  s(t  - - —  )  are  simply  delayed  versions  of  s(t)  and  c 

6_  •  x 

is  the  speed  of  sound  in  water.  Note  tnat  — —  is  the  time  delay  at 
point  x  relative  to  a  reference  point  at  the  origin,  with 


fT  12 

J-T/2 


s(t,x)  e  dt  ,  so  =  — ^ — 

'  -  ’  k  T 


(2.13) 


nTt 


k  =  1,  2,  ...,  N  , 
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the  vector  from  the  origin  to  x  defined  by  Defining  ^  -  t 
(2*13)  can  be  written  in  the  form 


S(Vx) 


6_  •  X_ 
C~ 

A  •  A 
~ 


&_  •  X 

Si*)  -p - 

Jt 


d* 


5  •  X 


6  •  X 


s(*>  *-=— 
JT 


Linder  assumption  4,  (2.14)  becomes 


where 


6_  •  _X 
—  j  U)  - 

S(«k,x)  =  S(u>k)  e 


“K1 

s(  t)  - dt  . 

Jr 


Now  define  sk(x) 


jv  x 

=  e  with 


a  uk 
vs  c 


cos  es;  then 


&_•  x 

- - —  » 

c 


(2.14) 


(2.15) 


(2.16) 
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SUk,x)  =  S(u)k)sk(x) 


(2.17) 


is  the  product  of  a  Gaussian  random  variable,  S(ui  ),  describing 
temporal  statistics  of  s(t)  and  a  deterministic  function,  s,(x), 

K 

representing  the  spatial  aspects  of  s(t,x)  at  the  kth  frequency. 

★ 

Since  E[Slu>k)S(u>  )  j  =  0  under  assumption  4,  the  space-time  covariance 
function  for  the  signal  is 


Rs(u,j<’ui  ;x>y)  =  e 


S(u>k,x)S(u>j  ,y ) 


=  l  E 


S(u>k) 


O  , 


jv$(x-y),  if  k  =  1 


if  k  £  1  • 


(2.18) 


Again  under  assumption  4, 


the  power  spectrum  of  the 


radiated  signal  process.  Fourier  transforming  (2.18)  witn  respect  to 
the  variable  a  =  x  -  y,  one  oDtains  the  wavenumoer  spectrum  for  the 
signal : 


S(i 


i  A 
,v)  = 


n  )e 


-jvQ 


dft 


do 


=  2tt  Sk  6(v  -  v$)  , 

(2.19) 


where  6(v)  is  a  Dirac  delta  function.  Tne  f requency-wavenumber 
spectrum  for  the  signal,  at  frequency  &>,  ,  is  thus  modeled  as  a  delta 

K 

function  at  v  =  v  with  known  amplitude  2n  S;  . 

S  K 
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Equation  (2.15)  characterizes  the  signal  at  x.  The  equivalent 
for  a  discrete  set  of  sensors  at  (x^,  x^,  . ..,  x^)  is 


/•T/2  .  . 

s(u>.  )  =  I  s(t)  e  dt  =  S(u),  ) 

At/2  ~  ir~ 


-Juk  5.  • 


"Juk 


,  (2.20) 


which  is  a  vector  of  essential  spatial  or  geometric  information  scaled 
by  the  random  variaDle  S(u>  ).  S(u)  )  has  the  covariance  matrix 

l\  K 


W  S  Sk  pk 


S  V  V 
K  \  -U<. 


(2.21) 


where  the  "steering  vector"  V,  is  defined  by 


A  • 

~~ c 


(2.22) 


i*in 
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Note  that  is  a  rank  one  matrix  with  a  single  nonzero 
eigenvalue: 

*1  =  trace  Pk  =  M  ,  (2.23) 


xi  =  o  i  1 


(2.24) 


2.4  Noise  Model 

We  assume  the  following  properties  for  the  noise  process  n(t,x)  = 
nc(t,x)  +  w(t,x),  where  nc(t,x)  represents  the  flow  component  and 
w(t,x)  represents  a  spatially  uncorrelated  noise  component.  These 
assumptions  are  largely,  but  not  entirely,  equivalent  to  our  earlier 
assumptions  concerning  the  signal. 

1.  The  noise  components  are  assumed  to  be  statistically 
independent  of  each  otner  and  of  tne  signal. 

2.  n(t,x)  is  a  sample  function  from  a  zero  mean  Gaussian  process 
characterized  by  known  second  moment  properties,  space 
homogeneity,  and  time  stacionari ty . 

3.  The  observation  time  T  is  large  compared  with  the  correlation 
time  of  tne  total  noise  process  n(t,x)  and  the  time  required 
for  nc(t,x)  to  travel  its  correlation  distance. 
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The  noise  spatial  covariance  function,  at  a  pair  of  frequencies,  is 
given  by 


where 


3 

V 

3 

C T 

X 

1 

<*< 

11 1> 

m 

i - 1 

rT/2  i 

rT/2 

dt 

da  E 

■'-T/2  J 

-T/2  L 

n(t,x)n*(o,y)  e 


>y)J 


-j'-»kt+jw^a 


T/2  yT/2 


dt  f 

T/2  •'-T/2 


da  R  (t-a;x-y)e 


-jw^x+ju^a 


(2. 


/•T/2  _•  . 

N  ( u>,, ,  x )  =  I  n(t,x)e  dt 

J-T/2  Jf 


(2 


and 


Rn(t-o;x-y)  =  E 


n(t,x)  n  (a,y ) 
Thus,  wnen  conditions  2  and  3  are  satisfied. 


(2 


q(u>K,wj  ;x-y) 


-ju>,  ( t-a ) 

da  Rn(t-a;x-y)e 


=  qk(x-y). 


if  k 


if 


25) 

.26) 

.27) 


K  *  1  . 
2.28) 
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Observe  that  the  noise  covariance  function  is  homogeneous  and  thus 
depends  only  on  the  separation  between  spatial  observations. 

The  power  spectrum  for  the  noise  at  a  point  is  given  by 

lim  q.  ( o )  =  litn 
T»co  T»a> 


T 12 

IS 

-T/2 


dt  do  R  (t  -  a;o)e 


-Juj^t-a) 


=  N, 


(2.29; 


The  wavenumber-frequency  power  density  for  tne  noise,  at  frequency 
u>k,  is  obtained  by  Fourier  transforming  the  stationary  covariance 
function  in  (2.28)  with  respect  to  n  =  x  -  y,  yielding 


w  _jQv  /-oo  A 12  A 12 

Q„  (v)  =  /q.  (a)e  d  £1  =  y  I  dal  dt  I  do  R  (t-o;a)e 

'  Jco  JAI2  J-TI2  " 


-j:ok  (t-o  )-jv£l 


K  =  1 ,  2,  ...,  N  . 


(2.30) 


The  physical  observations  discussed  in  Chapter  1  indicate  that 
the  flow  noise  component  of  Q,  (v)  has  tne  form  of  Fig.  2.3.  For 
analytical  purposes,  we  model  the  flow  noise  component  of  Q.  (v)  by 

K 

the  Butterworth  function: 


F  (  )  23|<  < 

Fk(v)  =  - p - 7 

(v  -  v J  +  a„w 


(2.31) 


K  =  1,  2, 
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Thus,  v,  characterizes  the  center  wavenumber  and  3.  tne  oandwidth  of  the 

K  K 

spatial  spectrum  with 


and 


a  ,  a  constant. 


(2.32) 


convective  speed 


(2.33) 


1^  is  the  flow  noise  spectral  level  at  frequency  The  flow 

noise  component  of  qk(n)  is  the  inverse  Fourier  transform  of  (2.32): 


Pk(Q)  *  *K 


-B<|0|+jvkQ 


(2.34) 


Equation  (2.31)  is  thus  equivalent  to  the  Corcos  [27],  [28] 
empirical  model  for  a  turbulent  flow  pressure  field,  which  is  a 
plausible  model  for  tne  T8L  noise  observed  by  a  towed  array. 
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If,  in  addition  to  the  flow  noise,  there  is  a  spatially  uncorrelated 
noise  component  of  spectral  level  W  . ,  then 

K 


q,  ( x  .  -  x  • )  =  I  .  e 
Hkv  i  j'  k 


-Sk|xrx^j 


w.  6..  , 

N  1  J 


(2.35) 


1,  2,  ...,  N  , 


where  6 . .  is  the  KronecKer  delta  (6)  function.  The  Kronecker  delta 
model  for  the  uncorrelated  noise  process  yields  finite  power  density  at 
eacn  oDservation  point.  Thus,  noise  covariance  matrices  have  finite 
traces . 

Observe  that  the  postulated  noise  model  in  (2.31)  allows  us  to 
adjust  all  of  the  typically  important  variaoles  (center  wavenumber, 
oandwidth,  and  power  level)  discussed  in  Chapter  1.  Because  (2.31)  is 
rational,  solutions  to  the  integral  equations  arising  in  the  binary 
detection  problem  are  obtained  in  a  straightforward  manner. 

The  uncorrelated  noise  makes  no  contribution  to  the  wavenumber 
spectrum  with  continuous  observations.  To  show  this,  we  rewrite  (2.35) 
in  the  form 


-s,  a 


qK(a)  = 


+jvkn 


wKf(o) 


(2.36) 
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where 


i,  if  «  =  o 


f(n)  = 


(2.37) 


0,  if  fi  o 


Fourier  transforming  (2.36),  we  obtain 


Ike  kl  1  J'  k  (2.38) 


For  finite  W.  ,  the  contribution  from  the  second  term  is  zero. 

As  long  as  one  is  dealing  with  continuous  observations,  only  the  flow 
noise  component  of  (2.35)  need  be  considered.  We  must  re-examine  this 
issue  wnen  arrays  of  discrete  sensors  are  discussed.  We  shall,  however, 
argue  for  ignoring  uncorrelated  noise  in  this  case,  as  well. 

When  there  are  M  sensors,  the  noise  wavesnapes  are  described  by  the 
M  vector  £(t)  or  the  correspondi ng  column  vector  of  Fourier  coefficients 


N(^k,x1) 

N(cok,x2) 


dt 


k  =  1,  2 


■  ■  •  » 


N  . 


(2.39) 
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Since  the  noise  processes  are  jointly  Gaussian,  the  vector  process  is 
described  completely  by  the  MN  x  iv1N  covariance  matrix 


(2.40) 


k  =  1,  2,  N  , 


consisting,  for  large  Tw,  of  Dlocks  of  M  x  M  matrices  of  the  form  in 

(2.12). 

2.5  Likelihood  Ratio. 

2.5.1  Realization  for  Discrete  Observations 

Well-known  analytical  results  [19]  assert  that  the  best  detector 
forms  the  likelinood  ratio 


(2.41) 


and  compares  it  to  a  threshold.  Pvr_/H^)  is  the  conditional  probability 


density  of  the  received  data  vector,  r,  in  (2.5)  under  hypothesis  H^. 


The  general  form  of  the  log  likelihood  ratio  for  Gaussian  signals 
and  noise  is  given  by  [19,  pp.  296—1 1 6 j .  For  the  complex  data  vector  in 
(2.5),  the  log  likelihood  ratio  becomes  [30] 
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=  In 

det 

i — H 

+ 

* 

+  r 

[k  -1  -  ( 

'k  +  K  )  -1! 

n\  n  sy 

L n  \ 

V n  s/  J 

£  > 
(2.42) 


where  the  covariance  matrix  under  the  noise-only  hypothesis,  from 
(2.40),  is  assumed  to  be  nonsingular  and 


with 


K 


n 


r  r*/h|1  , 


(2.43) 


K  =  E 


r  - 

£  £ 


/signal  only 


]■ 


(2.44) 


Making  use  of  the  dIock  diagonal  structure  of  the  covariance 
matrix  in  (2.12),  an  equivalent  test  statistic  is 


z  =  In 


[L<I>]  -  'nj^det(Kn(Kn  *  S>*1] 


£ 

k=l 


4 


<W*  -  <SA  *  W 


-1 


^k  ’ 


(2.45) 


where 


Q 


A 


k  " 


W'Nk 


(2.46) 
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and  is  defined  by  (2. PI), 


Appealing  to  Bartlett's  identity  [31], 


<\\  *  W1  “  lNkQK)_1 


2-1  *  -1 
N,  Q,  1  V,  V,  Q,  1 

k  y k  -k  -k  xk 

S 

k  *  -1 

1  +  — -  V  0  V 

1  N,  \  \  \ 


(2.47) 


Inserting  this  into  (2.45),  we  obtain 


z 


E 


S  /  N  ,/ 
k  k 


1  +  Ik  , 

N  ko 

k 


(2.48) 


where 


G 


G  (w,  )  =  V  Q  V 

o'  k'  -k  sk  -k 


(2.49) 


is  the  optimal  array  gain. 
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Significant  detection  problems  exist  primarily  at  low 
signal-to-noise  ratios.  With  SK / <<  1,  one  obtains  the 
small  signal  likelihood  ratio  test: 


z 


s 


A 


(2.50) 


If  Qk  and  the  spectral  functions  are  approximately  constant  over  a 
frequency  interval  of  length  (2i r/T),  the  k  sum  can  be  approximated  by 
an  integral: 


z 

s 


1_ 

Aw 


4 


2 

Aoj  > 


n/S  ( to ) 


V  (u) )  Q  ( u ) 


-1 


r(w) 


2 

doo  , 


(2.51) 

2tt 

with  Aw  =  — —  »  du  and  Nao>  >  w,  the  upper  end  of  the  processed 
frequency  band. 


An  implementation  of  the  small  signal  test  is  shown  in  Fig.  2.4, 
where  we  have  appealed  to  Parseval's  theorem  [21]  to  replace  the 
frequency  sum  in  (2.51)  by  a  time  sum. 
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2.5.2  Space-Continuous  Observations 

We  now  consider  a  heuristic,  but  direct,  argument  to  derive  the 
likelihood  ratio  test  for  space  continuous  observation  in  [-L/2,  L/2], 
The  more  rigorous  argument  is  made  in  [32]. 

With  the  definition 

f-k‘  'AhW1  •  (2-52) 


the  kth 


component  of  (2.50) 


can  be  expressed  as 


(2.53) 


wnere  the  notation  <•  ,  •>  denotes  an  inner  product  in  the  complex  M 
space,  C  .  Let  f^,  v^,  and  r.  .  be  the  jth  components  of 
f^,  V^,  and  respectively.  Now  define 


’kj 


=  1,  2, 


(2.54) 


39 


TO  7997 


where  a  is  the  uniform  spacing  oetween  oDservation  points.  TaKing  the 
transpose  conjugate  of  (2.52)  and  then  premultiplying  by  NkQk,  we 
have 

M 

Y''  N,  Q.  ■  g,,A  =  >Js, v  .  (2.55 

kxij  akj  < 

j=l 

and 


where 

i  =  1 ,  2 ,  . . . ,  M , 
k  s  1,  2,  . . . ,  N,  and 

M  =  x  +  1  , 

A 

with  Q.  .  the  ith  jth  element  of  Qk- 

For  a  »  o  (arbitrarily  tight  sensor  spacing  on  [-L/2,  L/2]),  with 
generalized  function  g  and  smootnly  varying  functions  q  and  r,  the 
sums  in  (2.55)  and  (2.56)  can  oe  approximated  by  integrals  on 
C— L/2,  L/2]  yielding 


E 

j-i 


kj  kj 


(2.56) 


qk(K,y)gK(y)dy  = 


-L/2  <  x  <  L/2 


(2.57) 
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and 


z 


s 


K 


2 


9k(x),  rk(x 


2 


(2.58) 


-L/2 


where 


\Qu  -  \<vyj> »  ^kCx’y>  • 


(2.59) 


k  *  1,  2,  N  . 


Equation  (2.58)  shows  that  the  likelihood  ratio  test  with  continuous 
observati ons ,  at  frequency  is  tne  magnitude-squared  inner  product 
between  the  data  r  (x)  and  the  function  g  (x).  To  construct  the 

K  K 


optimal  processor  and  thus  the  sufficient  statistic  z*  ,  we  need  only 


determine  g  (x)  from  (2.57)  for  our  Butterwortn  flow  noise  model  from 

r\ 

(2.34).  Since  we  consider  colored  noise  only,  tne  relevant  integral 
equations  are  Fredholm  integral  equations  of  the  first  kind. 

Under  the  normal  conditions  for  Riemann  integrabi 1 i ty,  the  sum  over 
k  implied  in  (2.58)  can  be  approximated  by  an  integral  for  sufficiently 
large  T: 


z 


s 


(2.60) 


c.  Tf 

where  A  u  =  — j  ■>  dw  and  N  Au  >  w  as  T  ■>  oo  . 
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The  more  rigorous  treatment  of  the  transition  trom  discrete  to 
continuous  observations  in  l 32 J  snows  that  the  solution  of  (2.55) 
converges  formally  to  the  solution  of  (2.57)  if,  and  only  if,  the 
likelihood  ratio  test  is  nonsingular,  a  condition  met  in  all  problems  of 
practical  interest. 

2.6  Likelihood  Ratio  Test  Performance 

2.6.1  Detection  Index  for  General  Detector 

When  the  input  covers  a  wide  oand  of  frequencies,  the  familiar 
detector  structure  in  Fig.  2.4  is  often  characterized  by  a  detection 
index,  d2,  defined  in  terms  of  the  output  quantity  z  by  the  equation 


M 

H 

Li 

K1 

2 

Var(z/H 

(2.61) 


2 

Nuttall  and  Hyde  [5]  have  obtained  a  general  expression  for  d  in  terms 
of  an  arbitrary  matrix  filter  function  K  : 


? 

d1"  = 


V. 

-K 


★ 

h  4 


(2.62) 
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The  narrowband  performance  of  the  general  detector  instrumentation  is 
widely  characterized  by  the  array  gain  Gk  at  frequency  .  In 
physical  terms,  G^  is  the  signal -to-noise  ratio  at  the  matrix  filter 
output  (Point  B  in  Fig.  2.4)  divided  by  the  signal-to-noise  ratio  at 
each  sensor.  Analytically,  it  takes  the  form  [33] 


Gk  = 


*  * 

H,  V  V,  H. 
-k  -k  -k  -k 


(2.63) 


-k 


gk  ** 


v  is  the  steering  vector  of  the  signal  from  (2.23).  Again, 

K 

★ 

is  the  normalized  spatial  covariance  matrix  of  the  received 


signal . 


If  Qk  is  positive  definite  and  is  the  largest  eigenvalue  of 


-1 


Qk  \  ,  a  well-known  theorem  of  linear  algebra  asserts  that 


<  —  X1  ’ 


(2.64) 


with  equality  if,  and  only  if,  H,  is  an  eigenvector,  flL , 

K  i. 

associated  with  x^.  Since  has  rank  1,  there  is  only 

one  nonzero  eigenvalue  and  0,  is  unique  to  within  a  multiplicative 
constant,  CK.  The  fact  that  multiplication  of  by  a  scalar 
constant,  C^,  does  not  alter  G^  is  immediately  obvious  from  (2.63). 
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\ 


Since 


Qk  1  Vk  VK*  has  rank  1  , 


=  trace  (Q"1  Vk  Vk*)  =  trace  Q'1  yk)  =  yk*  Qk_1  Vk  =  Gko 


(2.65) 


This  is  the  familiar  expression  for  the  optimal  array  gain  GkQ, 
first  encountered  in  (2.49).  The  eigenvector  0 ^  must  satisfy 


Qk_1  -k*  ^1  =  -k*  Qk  1  -k  h 


-1 


This  equation  is  satisfied  oy  0^  =  Qk  Vk  so  that 


Hklopt  =  ^  -k  ' 


un 


ique  to  within  a  multiplicative  constant,  C 


(2.66) 


(2.67) 


Instead  of  using  (2.67),  conventional  array  processing  procedures 
simply  align  the  signal  components  of  the  received  waveshape,  which 
amounts  to  the  choice 


^k|con  — k 


(2.68) 
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Since 


(2.63)  yields  the  conventional  array  gain 


kc 


v  *  q.  v, 

-k  wk  -k 


(2.69) 


All  of  these  results  are  familiar  from  the  literature  [33],  [34]. 
Because  of  the  invariance  of  6,  to  scale  changes  in  H  ,  one  does 
not  require  equality  of  (2.68)  and  (2.67)  to  achieve  equality  of 
conventional  and  optimal  array  gains.  Instead,  one  needs  only 


V1  \  -  ck  h 


(2.70) 


In  other  words,  the  steering  vector  must  be  an  eigenvector  of 
^  and  hence  of  Q^. 

If  the  sensor  separation  is  much  larger  than  the  noise 
correlation  distance  (1  «  &  a),  q,  is  an  identity  matrix  and 
(2.70)  is  satisfied  with  C,  =  M.  If  g.  a  <  1  but  s  L  >>  1,  it 

K  K  K 
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follows  that  M  >>  1.  Since  Q.  is  a  Toeplitz  matrix,  the  left  side 
of  (2.70)  has  elements 


Q,  1  V, 

Hk  —t 


Ql-n  6 


j2irpAvc  j2irlAvc 


1-M 


‘l-p 


Y  q. 


-j  2ttrAv 


1  p=l 


r=1-l 


s  , 


(2.71) 


where  r  =  1  -  p.  Qf  differs  from  zero  only  for  r 


3,,  A 


Hence,  the  r  sum  is  independent  of  1,  except  for  edge  effects  near 
1  =  1  and  M.  Except  for  tnese  edge  effects,  (2.70)  is  therefore 
satisfied  with 


C 


k 


-j2itrAvs 

e 


(2.72) 


Tnis  somewhat  informal  argument  suggests  tnat  the  condition  >>  1 
is  sufficient  to  ensure  at  least  approximate  equality  of  the  optimal 
and  conventional  array  gains  of  uniformly  spaced  linear  arrays. 
Precise,  but  less  general,  calculations  presented  later  will  confirm 
this  conclusion. 
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It  is  important  to  observe  that  a  frequency-dependent  rescaling  of 
by  a  scalar  will  not,  in  general,  leave  d^  unchanged. 

Thus,  different  instrumentations  may  have  the  same  array  gain  but 
dramatically  different  detection  indices. 

2 

According  to  a  well-known  result  in  detection  theory  [33],  d 
is  maximized  by  the  choice 


a*  -  if  V1  \  •  <2-73> 

K 

Thus,  in  addition  to  the  spatial  operation  (2.67),  the  best  processor 
uses  the  "Eckart  filter"  (see  Fig.  2.4),  with  transfer  function 
/ Nk *  Since  comillon  to  al  I  cnannels,  it  need  be  built  only 

once  and  is  inserted  just  ahead  of  the  squaring  operation  in  Fig.  2.4. 

/Jith  the  given  by  (2.73),  the  optimum  detection  index  becomes 


(2.74) 
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For  large  T,  one  can  make  the  usual  integral  approximation 


d_  = 


:T_fW 

”  2lIJ 

Jo 


V*U)Q(u>)  1  V(u»] 
N(u>)  ‘ 


du) 


/•w 

T_ 

2  TT 

L 

S(u>)  p  (  \ 
JR  V“> 

do>  , 


(2.75) 


where  again  G  (<d)  is  the  optimal  array  gain  at  frequency  to. 

Using  the  same  Eckart  filter  with  the  conventional  beamformer  in 

(2.68), 


(2.76) 


Equation  (2.62)  becomes 


(2.77) 


G  is  the  conventional  array  gain  given  by  (2.69). 
kc 
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The  sum  in  (2.74)  is  dominated  by  frequencies  where  the  optimal 

array  gain  G^q  (and  hence  the  postbeamforming  signal-to-noi se  ratio) 

is  large.  In  contrast,  the  denominator  sum  in  (2.77)  is  dominated  by 

frequencies  where  G  is  small.  Low  conventional  array  gain  at  some 

frequencies  will,  therefore,  dramatically  impair  the  performance  of 

the  conventional  detector  wnile  it  will  only  have  a  modest  impact  on 

the  optimal  detector.  This  statement  remains  true  even  when  G, 

ko 

G^c  over  the  entire  frequency  band  of  interest  (unless  the  two  array 
gains  are  frequency  independent).  Tne  only  possible  explanation  is 
that  the  simple  Eckart  filter  is  well  matched  to  the  optimal  processor 
but  poorly  matched  to  the  conventional  oeamformer.  The  ideal 
postbeamformer  filter  (PBF)  should  be  chosen  according  to  its  input 
signal  and  noise  signal  and  noise  spectra,  not  the  element  signal  and 
noise  spectra.  The  latter  choice  employs  'n/S^/N^  as  indicated  in 
(2.76).  We  must,  therefore,  re-examine  the  issue  of  proper  scalar 
filcering  after  tne  initial  choice  of  the  spatial  filter  has  been  made. 

We  now  modify  the  detector  structure  in  Fig.  2.4  by  introducing 
the  scalar  filter  C(w),  as  shown  in  Fig.  2.5. 

In  what  follows,  we  keep  JH ( u> )  temporarily  fixed  and  adjust  C(w) 
for  the  optimum  detection  index.  Equation  (2.62)  now  reads 
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k 

Ck 

n  *  * 

s,  H.  V,  Vu  H, 

K  -k  -k  -k  -k 

2 

£ 

k 

CJ 

4  N  2 

1  K 

2 

(2.78) 


has  simply  been  replaced  by  CK 


Theorem: 


Max  d  = 

C, 


£ 


V.  2 

N,  2  GK 

K 


(2.79) 


where  Gk  is  given  by  (2.63)  and  the  maximizing  C;<  satisfy 


C„ 


2  Sk 


N 


2 


: k 

2 

^  Qk4 

(2.80) 


Proof: 


Let  \  E  »k*  *k  -<*  «k 


Bk "  \  ^k  \  • 


(2.81) 

(2.82) 


50 


TR  7997 


Then  (2.78)  reads 


(2.83) 


Rearranging  the  numerator  and  using  the  Scnwarti  inequality,  one 
obtains 


with  equality  if,  and  only  if, 


(2.85) 


where  C  is  a  frequency  invariant  constant.  Substituting  (2.85)  into 
(2.83)  and  using  (2.81),  (2.82),  and  (2.63), 
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4*  h  h  4 


4 


(2.86) 


Since  the  frequency  invariant  constant  C  does  not  affect  the  detection 
index,  the  maximizing  |ck |2,  i.e.,  the  postbeamformer  filter  power 


transfer  function,  can  be  written  as 


C„ 


*  * 

4  4 


★  1 2 


(2.87) 


This  completes  the  proof. 

With  the  optimal  spatial  filter  *  V^,  (2.87) 

becomes 


(2.88) 


Thus,  we  confirm  that  the  simple  tckart  filter  is  indeed  optimal  when 

matrix  filtering  is  optimal.  On  the  other  hand,  for  the  conventional 

beamformer  H,  =  V,  ,  the  best  frequency  filter  power  transfer 
K  K 

function  is  given  by 
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r 

2  k 

M'- 

1 

K 

’  \2 

v/  Q  v, 

-k  K  — K 

2  "m2 

r  z 
,  2  bkc 

V 


(2.89) 


2 

The  frequency  invariant  constant  1/M  is,  of  course, 
immaterial.  The  important  feature  of  (2.89)  is  the  last  factor:  the 
best  frequency  filter  voltage  transfer  function  requires  the  weight 
in  addition  to  the  Eckart  filter.  With  this  weighting,  the 
dependence  of  the  detection  index  on  the  input  signal-to-noise 
spectral  ratio  and  array  gain  (2.79)  is  precisely  the  same  for  the 
optimal  and  conventional  detector  (or,  for  that  matter,  any  detector 
with  specif ied  spatial  filter  rt  ) .  Differences  in  performance  arise 
only  because  of  differences  between  the  optimal  and  suboptimal  array 
gains,  i.e.,  choices  of  H.  . 

The  above  ideas  are  easi ly  extended  to  space-continuous 
observations.  Using  (2.52),  we  have  from  (2.74) 


d 

o 


2 


K 


k 


* 

f,  V. 
— k  — k 


(2.90) 
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With  the  notation  for  gkj  in  (2.54),  the  continuous  formulation  of 
optimal  wideband  detection  performance  follows  easily: 


with  gk ( x )  »  g(u>,x),  >  S(u>),  and  >  w. 

Equation  (2.91)  says  that  optimal  wideband  detection  performance 
in  the  continuous  limit  is  specified  completely  oy  the  inner  product 

jv  X 

of  the  signal  model  e  and  the  function  g(u),x),  whicn  one  obtains 
directly  from  (2.57).  The  inner  product  in  (2.91)  is  over  [-L/2,  L/2]. 
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y(t,x)  =  n(t,x)  +  s(t,x) 


Fig.  2.1.  Aperture  With  Plane  Wave  Signal  and  Additive  Colored  Noise 
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Fig.  2.2.  Signal  Model 
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C  =  SPEED  OF  SOUND 


Vs  =  27rfk ,  c  cos  0, 


lie  +jMl 

n  =  \  -  y 


—  2 tv  Av  cos  Os 


Fig.  2.3.  Coherent  Wavefront  Signal  and  First  Order-Butterworth 
Wavenumber  Spectrum  at  kth  Narrowband  Frequency 
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M  SENSORS 


FILTER  OUTPUT 
(SCALAR) 


INTEGRATOR 


H(u>)  =  V  (w)  QM 


MATRIX  FILTER  FOR:  _ 

1)  MAX  S/NAT  A:  v  S-lllL  ECKART  FILTER 

N(a>) 

v  s(co) 

2)  BEST  FIT  TO  SIGNAL  WAVESHAPE  AT  A:  - - 

N(oj) 


Fig.  2  4.  Small  Signal  Likelihood  Ratio  Test 


M  SENSORS 


SCALAR  FILTER 


INTEGRATOR 


Fig.  2.5.  Detector  for  Random  Signal  and  Noise  With  Frequency  Scaling 
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CHAPTER  THREE 


CONTINUOUS  SPATIAL  OBSERVATION 


3.1  Optimal  Processor  Structure 

Tnis  chapter  deals  with  detection  in  a  flow  noise  environment, 
using  a  continuous  observation  interval  [-L/2,  1/2]. 

Van  Trees  [19,  p.  316]  describes  a  general  technique  of  solving 
(2.57)  for  processes  with  rational  spectra.  If  the  noise  spectrum  has 
the  form  of  (2.31),  his  procedure  yields 


D  C 

+  -|  6(x  +  L/2)  +  s(x  -  L/2) ,  — L / 2  <  x  <  L/2 

(3.1) 


where 


and 


A 


(3.2) 


V  ,  =  V,  -  V 

d  k  s 
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Substitution  into  (2.57)  easily  verifies  that  (3.1)  indeed  satisfies 
the  integral  equation. 

According  to  (2.58),  one  constructs  the  spatial  part  of  the 
likelihood  ratio  test  at  frequency  a  as  an  inner  product  between 

K. 

rk(x)  and  gk(x) : 


(x)rk 


( X )  dx 


*  * 
rK(-L/2)  +  \  r(+L/2)  + 


+ 


-jv$x 

dx  r.,(x)e 


(3.3) 


The  structure  is  snown  in  Fig.  3.1.  *le  see  that  it  is  divided  into 
two  parts:  a  spatial-matched  filter  ^conventional  beamformer),  which 
correlates  the  received  process,  r  (x),  with  a  signal  replica  (i.e., 
s  *(x))  and  two  end-point  sampling  sensors,  which  receive  the 

K 

pressure  field  at  tne  interval  end  points.  The  scalar  Vs”k/Fk(vs) 
represents  tne  relative  weighting  necessary  for  the  conventional 
processor  to  achieve  near  optimal  wideband  performance. 
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3.2  Narrowband  Performance  of  the  Optimal  Spatial  Processor 


Nothing  is  lost  in  our  understanding  of  narrowband  optimal 

? 

performance  if  we  consider  only  the  optimal  array  gain  in  d  from 
(2.91): 


A  \  f  -jV 

Gkn  =  ~F~  I  9kU)  e  dx,  k  =  1,  2,  N  (3.4) 

K0  Vs  J- L/2  K 


dq  jv$L/2 
2  6 


c0  -jvsL/2 
2  e 

(3.5) 


Substituting  for  bQ  and  cq  from  (3.2),  one  obtains 


V  ■  +  1 


,n+V  fi  +  v« 


(3.6) 


This  result  is  intuitively  appealing.  If  the  aperture  length  L 
approaches  zero,  then  the  optimal  instrumentation  becomes  effectively 
a  single  sensor  with  unity  array  gain: 


G  >  1 

KO 


(3.7) 
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On  the  other  hand,  if 


then 


1kL 


»  1 


(3.8) 


(3.9) 


which  provides  the  intuitively  satisfying  ooservation  that 


Total  noise  power 


ko 


W111 


Noise  power  in 


band  of  width  1/L  near  v. 


(3.10) 


must  exceed  unity  for  the  optimal  array  gain  to  exceed  unity.  In  the 
absence  of  wavenumber  structure  for  the  noise  in  in 

Fig.  2.3.,  i.e.,  «  1,  optimal  array  gain  significantly 

exceeds  unity  if,  and  only  if,  exceeds  2. 

In  the  absence  of  white  noise,  one  can  achieve  perfect 
detectability  (singular  detection)  [19,  p.  303]  with  optimal 
processing,  namely 
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G 


»  oo  , 


(3.11) 


KO 


if  the  spatial  bandwidth,  s  ,  is  exactly  zero  in  (3.9).  For  flow 
noise,  cannot  be  zero  and  therefore  this  problem  will  not  arise. 

The  result  in  (3.9)  is  very  clear.  Array  gain  with  optimal 
processing  is  large  if,  and  only  if, 


(3.12) 


For  — 7-  >>  1,  tnis  inequality  always  holds.  For  — <  1,  it  holds 

only  if  the  noise  wavenumber  spectrum  is  narrow  enough,  i.e.,  the 
fractional  bandwidth  is  small  enough. 


There  is  one  additional  point  worth  making.  If  v^L  <  1,  we 


shall  demonstrate  that  a  conventional  beamforiner  cannot  separate 
signal  from  noise  and  has  an  array  gain  near  unity.  The  condition 
<  1  with  large  enough  v^/6^  in  (3.6),  therefore,  appears  to 
offer  a  possibility  for  the  optimal  detector  to  outperform  the 
conventional  detector  in  this  case. 
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3.2.1  Wideband  Optimal  Detection  Index 

A  general  expression  for  the  wideband  optimal  detection  index  is 
obtained  directly  from  (2.90): 


(3.13) 


dQ2  is  simply  the  sum  of  weignted  optimal  array  gain  from  (3.6)  and  is 


therefore  most  heavily  influenced  oy  frequencies  where  the  product 


is  largest,  i.e.,  the  beamformer  output  signal-to-noise  ratio  is 
highest. 

Previously  discussed  properties  of  array  gain  obviously  carry 
over  directly  to  the  wideband  detection  index.  Improvement  over 
single  sensor  performance  is  clearly  substantial  if  ( 3 . 1 2 j  is 
satisfied  at  most  frequencies. 
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3.3  Narrowband  Performance  of  the  Conventional  Detector 


Tne  array  gain  for  the  conventional  processor  is  derived  in 
Appendix  A: 


1  + 


(vV) 

skl(\/V> 


~3kL 

e  (cos  v^L 


,  2vdBk  .  ,  x 

+  — - ^  sin * *  vdL^ 

vj  -0i, 


Gko  -  1 


-3kL  2v  ,S(/ 

e  (cos  v^L  +  — . . 2  sin  vdL) 

vd  “ek 


(3.14) 


The  near  optimality  property  of  conventional  processing  for 

e^L  »  1  follows  immediately.  It  is  also  clear  that  there  are  two 

distinct  mechanisms  of  array  gain  degradation  in  (3.14).  The 

"6kL 

sinusoidal  terms  in  (3.14)  vary  as  e  and  can  therefore  vanish 
quick  ly  when  b^L  >  1 . 


The  term 


2 

vd 


(3.15) 


55 


TR  7997 


varies  as  (s^L)  1  tor  vd  >  Bk  and  can  therefore  cause 


degradation  of  conventional  array  gain  at  b^L  values  appreciably 
less  than  unity.  Vie  show  later  that  conditions  s^L  <<  1  and 
V(j  <  Bk,  in  (3.14),  yield  Gkc  «  1. 

There  is  an  interesting  observation  to  be  made  from  (3.14).  The 

numerator  of  Gkc  is  essentially  equal  to  GkQ  if  \>dL(vd/23k)  >>  1. 

The  fact  that  GkQ  can  be  much  larger  than  Gkc  for  akL  <<  1 

indicates  that  the  end-point  sensors  must  achieve  noise  cancellation 

when  6.  L  <<  1  and  thus  are  responsible  for  improving  performance 
k 

over  the  performance  with  conventional  beamforming  alone. 

For  v  ,  »  S,  ,  the  noise  wavenumber  spectrum  in  |V»  Vi,] 
d  k  s  k. 

is  very  structured.  Equation  1 3 . 14 )  then  reduces  to  the  approximate 
form 


v  ,L  v, 
d  d 


2  B, 


(3.16) 


G, 


kc 


1  + 


If  also 
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so  that  the  optimal  array  gain  is  appreciaDly  larger  than  unity. 


This  ratio  clearly  exhibits  minima  near 


(3.17) 


v  .L  =  2nn,  n  =  1,  2, 


and  maxima  near 


\>dl  =  (2n  +  1 ) it ,  n  =  0,  1,  2,  - 

At  the  minima,  conventional  array  gain  in  (3.17)  is,  at  most, 

3  decibels  below  optimal  if  s^L  <<  1.  At  the  maxima  of  (3.17),  with 
<<  1,  conventional  array  gain  is  far  below  optimal: 


6kL 

2 


«  1 


(3.18) 


Fig.  3.2  shows  the  maxima  and  minima  of  (3.17)  plotted  as  a 
function  of  a,K L .  The  following  ooservations  are  plain.  There  are 
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no  more  than  3  decibels  of  advantage  with  optimal  processing  over 
conventional  processing  when  v^L  =  2iitt,  n  =  1,  2,  ....  Tne 
advantage  with  \>dL  =  (2n  +  l)ir,  n  =  0,  1,  2,  ...,  is  large  only  for 
8kL  <<  1  Dut  is  no  more  than  3  decibels  when  s^L  is  only  slightly 
larger  than  unity.  If  s^L  _>  2,  conventional  array  gain  is,  at  most, 
1.8  decibels  Delow  optimal.  The  practical  implication  is  that  the 
end-point  sampling  sensors  in  Fig.  3.1  provide  insignificant  array 
gain  improvement  when  the  continuous  aperture  L  exceeds  two  noise 
correlation  lengths.  Thus,  the  long  array  length  reduces  the 
detection  problem  to  near-white  noise  in  a  band  of  order  1/L  over 
which  the  response  of  the  conventional  beamformer  is  significant. 

This  result  is  intuitively  logical.  If  the  noise  is  near-white  over 
the  interval  where  the  beamformer  response  is  significant,  then 
regardless  of  what  processing  is  applied  at  the  beam  output,  it  should 
not  significantly  improve  performance. 

If  both  v^L  and  became  small,  one  would  expect 
conventional  array  gain  to  approach  unity  because,  then,  all  the  noise 
power  is  confined  to  the  main  lobe  of  the  beam  pattern.  To  establish 
this  fact  more  formally,  we  rewrite  (3.14)  and  expand  tne  nonalgebraic 
terms  into  an  appropriate  Taylor  series,  retaining  terms  to  second 
order  only: 
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26, 


K  +K2) 


2^V) 


,  J  2 

i  i_  , 
*  ) 


-i  V  +  3 
d 


2  2  v  ,2L2 

nr  Uvo  Vhi  --V)  - 


o  <  6  L  <<  1  , 

K. 

v  .L  <<  1  . 
u 


(3.19) 


Combining  terms  of  tne  same  order  in  L, 


yields 


2  ,2, 

2(v.j  -  b ^  J 

<»dz  * 


2^2  ~  \z) 


=  0 


(,3.20) 


26, 

\ 


23, 


(VH~  +  81,  )  (va  +  6,  ) 


=  0 


(3.21) 


L  , 


(3.22) 


and 


SKC  "  1  '  °  <  SkL  «  1  ' 


<<  1  , 


(3.23) 


1  i  m  G,  =  1  .  (3.24) 

K  L. 

L  >  o 
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In  the  limit  of  L  >  o  (fixed  sk  and  v^),  the  optimal  array 
gain  must  also  approach  unity,  a  fact  immediately  apparent  from 
(3.6).  However,  when  ek  is  small  enough  so  that 


v  ,L  v  . 
d  d  , 


the  optimal  array  gain  can  be  significant  even  though  the  conventional 
array  gain  has  fallen  close  to  unity.  Admittedly,  this  special  case 
may  have  only  academic  value  for  towed  array  applications  where  \>dL 
is  generally  mucn  larger  than  unity. 

These  observations  suggest  the  following  strategies  for  choosing 
aperture  lengtns  to  achieve  near  optimal  conventional  array  gain 
performance: 

1.  If  \>d  is  known,  choose  L  =  2 mi/vd  witn  n  as  large  as 
practically  possible.  Then,  conventional  performance  is  no  more  than 
3  decibels  below  optimal. 

2.  If  vd  is  not  known  exactly,  again  choose  L  =  2irn/\>d  with  n 

satisfying  v,/2mT  <  n.  But  this  is  equivalent  to  choosing 
Q  < 

L  >  1/s.  ,  and  hence  all  one  really  needs  is  an  array  longer  than  the 
noise  correlation  length.  Fig.  3.2  then  assures  one  of  a  worst  case 
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loss  of  <  3  decibels.  If  one  can  increase  the  array  length  to 
L  2  2/Bk,  the  worst  case  loss  is  reduced  to  1.8  decibels. 

3.3.1  Wideband  Conventional  Detection  Performance 

Array  gain  in  (2.63)  is  often  used  to  compare  performance  of 
alternate  array  processing  procedures.  It  specifies  postbeamformer 
signal-to-noise  ratio,  at  a  given  frequency,  in  terms  of  the  input 
signal-to-noise  ratio  and  is,  therefore,  independent  of  any  linear 
filter  connected  to  the  conventional  beamformer  output  (i.e., 
frequency  scaling  in  fig.  2.5). 

From  the  arguments  made  in  section  2.6,  we  Know  that  wideband 
detection  performance,  on  the  other  hand,  is  strongly  dependent  on 
such  a  linear  filter.  Hence,  two  detector  instrumentations  having 
identical  array  gains  may  have  very  different  wideband  detection 
indices.  We  also  know  that  oasic  Eckart  filtering  is  poorly  matched 
to  the  conventional  beamformer. 

Much  of  the  performance  loss,  as  expressed  by  (2.77),  incurred 
with  conventional  beamforming  followed  by  basic  Eckart  filtering 
(2.76),  can  be  recovered  by  the  simple  modification  of  the  basic 
Eckart  filter  (2.89).  The  expression  for  the  space-continuous 
conventional  detection  indices  with  oasic  and  modified  Eckart 
filtering  are  obtained  directly  from  (2.77),  (2.89),  and  (A. 4): 
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d 


c 


2 


simple 

Eckart 


d 


c 


2 


modified 

Eckart 


(3.25) 


(3.26) 


where  G  is  the  conventional  array  gain  on  continuous  [-L/2,  L/2] 

KC 

from  (3.14).  Again  the  variation  of  Gkc  in  (3.25)  witn  temporal 

frequency  causes  the  basic  Eckart  filter  to  perform  poorly  with 

conventional  beamforming.  Tnis  variation  is  illustrated  in  Fig.  3.3, 

using  the  Corcos  models  for  and  from  (2.32)  and  (2.33)  with 

v  =  0.  Variations  in  array  gain  with  frequency  are  directly 
s 

traceaole  to  the  fact  that  a  given  aperture  length  equals  an  even 
multiple  of  a  noise  naif  wavelength  at  some  frequencies  and  an  odd 
multiple  of  a  half  wavelength  at  other  frequencies. 

Fig.  3.4  snows  conventional  detection  indices  computed  from 

(3.16),  (3.25)  and  (3.26)  graphed  together  with  the  optimal  detection 

index  from  (3.6)  and  (3.13).  Parameter  values  are  =  a  =  0.01 

and  v  =  0.  Detection  indices  are  plotted  against  a  normalized 
s 

highest  signal  frequency,  wL/Uc.  Each  curve  is  obtained  by 
normalizing  with  respect  to  L  and  an  assumed  uniform  signal -to-noise 
ratio  in  the  signal  band. 
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Specif ical ly , 


Sk  S/I  (a  constant),  for  k  =  1,  2,  N 

k  o,  otherwise. 


(3.27) 


We  see  that  wideband  detection  performance  with  conventional 
beamforming  followed  by  modified  Eckart  filtering,  as  expected,  is 
everywhere  better  tnan  wideband  detection  performance  with 
conventional  beamforming  followed  by  basic  Eckart  filtering. 

As  wL/Uc  increases  to  2n,  one  encounters  the  first  conventional 
array  gain  peak  (see  Fig.  3.3).  Hence,  from  (3.26),  modified  Eckart 
filtering  oegins  to  recover  a  large  fraction  of  the  performance  loss 
incurred  with  oasic  Eckart  filtering. 

In  that  range  the  conventional  detector  curve,  with  the  modified 
Eckart  filter,  lies  a  fairly  constant  13  decibels  below  the  optimal 
detector  curve.  On  tne  other  nano,  the  performance  loss  of  the 
conventional  detector  with  simple  Eckart  filtering,  relative  to 
optimal,  ranges  from  38  decibels  at  wL/Uc  =  2-rr  to  48  decibels  at 
wL/Uc  =  1 2 n .  Modification  of  the  basic  Eckart  filter  thus  recovers 
a  suDstantial  fraction  of  the  performance  loss  even  when  ukL  <  1. 

Note  that  improvement  with  modified  Eckart  filtering  over  basic 
Ekcart  filtering  is  only  2  deciDels  when  wL/Uc  =  it. 
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One  expects  limited  benefit  from  the  modified  tcxart  filter  in  this, 
case  since  the  variation  in  array  gain  ^.Fig.  3.3)  is  relatively  small 
in  this  frequency  range.  In  this  extreme  low-pass  case,  use  of  the 
modified  tcxart  filter  obviously  becomes  less  important.  One  could, 
therefore,  consider  using  only  tne  simple  Ecxart  filter,  x/SK/IK, 
with  the  conventional  beamformer. 

One  could  raise  the  point  that  a  simple  high-pass  filter 

preceding  the  squaring  operation  in  Fig.  2.5  might  be  used  to 

eliminate  frequencies  where  array  gain  is  smallest,  i.e.,  the  interval 

[0,ir]  in  Fig.  3.3.  Since  ^3.25)  is  lowered  drastically  by  low  values 

of  G  in  even  small  frequency  bands,  elimination  of  such  bands 
KC 

might  secure  much  of  the  benefit  availaole  from  the  more  complex 
modified  Ecxart  filter.  This  suboptimal  approach  has  interesting 
practical  implications  and  is  pursued  in  Chapter  5.  The  high-pass 
filters  are  relatively  easy  to  construct  and,  as  we  shall  see,  are 
competitive  with  the  best  postbeamformer  filters  in  (2.89). 

When  b  L  »  1,  G  =  b  and  the  only  issue  then  is  the 

X  KO  XC 

benefit  to  be  gained  by  using  the  modified  Ecxart  filter.  From  (3.6) 
and  (3.14),  we  have 
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Vk  >>  VS 


(3.28) 


Using  spectral  functions  satisfying 


S  (to) 


j  (a  constant),  for  <_  w  _< 


0,  elsewnere 


(3.29) 


and  assuming  T(u>^  -  u, )  »  2n  so  tnat  sums  in  (3.25)  and  (3.26) 
can  be  approximated  by  integrals,  one  obtains 
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modified 

Eckart 
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simple 

Eckart 


2 tt  j2  (u2  ■  “r 
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2.2  2 

a  L  u) 
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- 2 - fv“ 2“  wi’“iu,2 


(3.31) 
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The  ratio  of  the  two  detection  indices  is,  therefore, 


modified 

Eckart 


simple 

Eckart 


(3.32) 


The  divergence  of  (3.32)  for  >  o  is  due  to  the  fact  that 

r  »  o  as  u),  >0  [see  (2.32)  and  (2.33)]  so  that  the  large  &.  L 
k  k  k 

assumption  is  violated  at  the  lower  end  of  the  spectrum. 

Significant  improvements  in  conventional  detection  are  obviously 
available  for  wide  bandwidth  signals.  For  instance,  if 


the  gain  in  (3.32)  is  approximately  5  deciDels,  Dut  if  uj u>^  =  50, 
the  gain  is  approximately  12  decibels. 

According  to  (3.28),  the  required  postbeamformer  filter  power 
transfer  function  is  characteri zed  by 
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2  Sk  2 
k  ~  ,  2  “k 

K 


(3.33) 


Thus,  with  b^,L  >>  1,  the  only  addition  to  the  simple  Eckart  filter 
is  a  voltage  gain  increasing  linearly  with  frequency.  The 
modification  is  independent  of  Uc  and,  therefore,  of  specific 
operating  conditions. 

3.4  Summary 

This  chapter  has  compared  tne  performance  of  conventional  and 
optimal  detectors,  using  space-continuous  observations.  The 
narrowband  comparison  is  summarized  in  Fig.  3.5.  Results  with  the 
Butterworth  flow  noise  model  indicate  that  one  can  achieve  near 
optimal  narrowband  performance  with  only  conventional  processing  in 
many  cases  of  practical  interest,  especially  if  ekL  >>  1.  An 
important  observation  is  that  one  achieves  this  performance  with  an 
instrumentation  that  is  relatively  simple  compared  with  the  fully 
optimal  construction.  The  ekL  >>  1  condition  is  certainly 
reasonable  for  most  practical  applications  where  overall  array  lengths 
usually  exceed  the  flow  noise  correlation  length. 

Over  such  a  space-continuous  length,  even  conventional  steering 
(beamforming)  is  impractical.  The  discussion  here  thus  serves  two 
purposes.  First,  it  sets  absolute  bounds  on  performance  that  can  be 


77 


TR  7997 


approacned  by  arrays  of  tightly  packed  discrete  sensors.  Secondly,  it 
leads  into  a  later  discussion  of  continuous  sensors. 

Wideband  detection  performance  is  strongly  dependent  on  linear 
filtering  applied  to  a  conventional  beainformer  output.  Two  detector 
instrumentations  having  identical  array  gains  can  therefore  have 
dramatically  different  detection  indices.  Basic  Eckart  filtering  is 
poorly  matched  to  the  conventional  beainformer,  especially  in  the 

<  1  regime.  A  simple  modification  of  the  Eckart  filter  can  very 
substantially  enhance  the  performance  of  the  detector  using 
conventional  beamforming,  making  it  more  competitive  with  the  optimal 
detector.  Major  gains  are  possible  with  modified  Eckart  filtering  and 
continuous  observations  when  the  noise  has  sufficient  spatial 
structure  or  a  reasonably  wide  band  of  frequencies  is  processed. 

The  modified  Ec<art  filter  instrumentation  depends  on 
conventional  array  gain,  which  is  not  usually  known  a  priori.  In  the 
large  regime,  the  instrumentation  reduces  to  a  relatively  simple 
pre-emphasis  of  high  signal  frequencies.  In  general,  array  gain  is 
not  likely  to  be  known  in  advance.  Hence,  an  adaptive  modified  EcKart 
filter  is  suggested.  However,  one  needs  only  a  single  adaptive  filter 
at  the  beam  output.  For  flow  noise  encountered  by  towed  arrays,  the 
primary  unknown  parameter  would  probably  be  the  convective  noise 
speed.  Adaptation  could,  therefore,  be  a  relatively  simple  matter. 
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In  Chapter  4,  we  revisit  the  binary  detection  problem  in  the 
context  of  the  sampled  aperture  and  examine  detection  issues 
surrounding  spatial  sampling,  such  as  the  attendant  detection  loss  due 
to  spatial  al i asing. 
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Fig.  3.1.  Optimal  Spatial  Structure  for  First  Order  Butterworth 
Wavenumber  Spectrum  and  Plane  Wave  Signal 
at  kth  Frequency 


80 


TR  7997 


Fig.  3.2.  Ratio  of  Optimal  to  Conventional  Array  Gain  at 
kth  Frequency  (vd/  #»1)  With  vd  =  yw  -  vs 
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Fig.  3.3.  Optimal  and  Conventional  Array  Gain  Versus  a>kL/Uc, 
a  =  /3k  vk  =  0.01,  ^kL^I  >  Vi  0 
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Fig.  3.4.  Wideband  Detection  Indices  Versus  wL/Uc  for  a  =  /3k/vk  =  0.01 
Normalized  With  Respect  to  SNR  (S/I  a  Constant), vs  =  (ws/c)cos  0S  =  0 
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Fig.  3.5.  Summary  of  Optimal  (Gk„)  and  Conventional(Gkc) 

Array  Gain  Performance 
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CHAPTER  FOUR 


DETECTION  PERFORMANCE  WITH  A  UNIFORMLY  SAMPLED  SPATIAL  APERTURE 

4.1  Introduction 

Continuing  tne  formal  analysis  from  Chapter  3,  we  specify  the 
optimal  detector  realization  and  corresponding  optimal  performance 
with  discrete  observations  regularly  spaced  in  [-L/2,  L/2]. 

4.2  Effects  of  Spatial  Aliasing 

Following  standard  procedures  of  sampled  data  analysis  [35,  pp. 
117-120;  36;  37],  the  spectrum  of  the  sampled  noise  process  can  be 
written  as 


00 


(4.1) 


J=-oo 


This  aliased  noise  wavenumber  spectrum  arises  from  taking  values  of 
the  homogeneous  Butterworth  covariance  function  in  (2.35)  at  discrete, 
but  uniform,  intervals  in  (-cd,cd).  The  wavenumber  density  for  the 
sampled  process  is  a  periodic  function  in  v,  with  period  2tt/a  obtained 
by  summing  translates  of  F^(v).  Since  the  input  noise  spectrum  is 
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not  zero  for  jvj  >  n / a  but  covers  an  infinite  number  of  wavenumbers, 
the  spectrum  of  the  sampled  process  will  necessarily  involve 
overlapping  of  individual  terms  in  the  summation.  In  all  cases,  the 
overlapping  results  in  increased  noise  spectral  levels  and  in  many 
cases  smoothes  the  noise  spectral  features.  One  obtains  the  mean 
square  value  for  the  sampled  process  by  integrating  the  sampled 
spectrum  over 


jl 

A  A  _ 

or,  more  generally,  over  any  period  of  length  2tt/a.  This  then  results 
in  finite  power  even  when  there  are  delta  functions  in  the  sampled 
spectrum. 

Notice  that  the  spectrum  for  the  sampled  process  equals  the 
original  spectrum  if,  and  only  if,  a  is  made  extremely  small  so  that 
discrete  observations  pass  in  the  limit  to  continuous  observations. 

For  our  Butterworth  wavenumber  model  (in  Fig.  2.3),  the  exact 
level  of  the  aliased  flow  noise  at  is  found  by  simply 
substituting  (2.31)  into  (4.1).  One  can  then  express  the  ratio  of  the 
sampled  flow  noise  spectrum,  at  the  signal  wavenumber,  to  the  flow 
noise  spectrum  with  continuous  observations  as 
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(4.2) 


If  the  fractional  spatial  Dandwidtn  is  assumed  small,  i.e., 
e,K/\>d  <<  1,  (4.2)  oecomes 
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(4.3) 


The  series  in  (4.3)  converges,  a  fact  easily  established  by  the  Limit 
Comparison  Test  of  elementary  calculus. 

It  follows  from  (4.3)  that 


2 
V 

-^pr  =  1  for  some  j  .  (4.4) 

Equation  (4.3)  may  be  used  to  compute  the  increase  in  noise  level  at 
v  caused  by  regular  sampling.  From  (4.4),  this  increase  is  clearly 
very  large  if  the  sampling  interval  a  is  an  integer  multiple  of 
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which  approximates  the  noise  wavelength 


2tt 


x,  =  — ,  when  v,  >>  v 
x  v.  x  s 

\ 


With  a  =  2-rr/v^,  the  nigh  wavenumber  noise  folds  down  (aliased)  onto 
the  signal,  maxing  it  impossible  to  separate  the  two  by  spatial 
processing.  Analogous  with  sampling  in  time,  noise  reconstruction 
from  the  sampled  values  would  loox  lixe  low  wavenumoer  acoustic 
components.  Thus  signal-to-noise  ratio  is  drastically  reduced 
relative  to  continuous  observations.  At  frequencies  for  which  A  is  a 
multiple  of  2ir/v  ,  noise  adds  coherently  between  sensors.  Noise 

K 

components  at  these  frequencies  can  easily  dominate  beam  output  noise 
power. 


By  contrast  when  a  =  (4.3)  yields 


F 


xa 


r„(v 

N 


s 


T 


=  2.4 


(4.5) 


which  is  an  increase  in  noise  level  at  vs  of  approximately  3.8 
decibels  relative  to  continuous  observations.  We  shall  show  that  this 
increased  noise  level  is  al 1  that  happens  in  the  transition  from 
continuous  to  discrete  observations.  Hence,  the  best  array  gain,  with 
discrete  observations,  at  this  spacing,  is  only  a  factor  of  2.4  [3.8 
decibels]  below  the  best  array  gain  with  continuous  observations. 
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The  conclusions  inferred  from  \,4.3)  are  of  fundamental  importance 
in  the  towed  array  application  where  v,  is  usually  much  larger  than 
vs-  One  avoids  the  aliasing  proolem,  in  this  case,  with  an 
extremely  tight  sensor  sampling  interval.  Tne  aliasing  problem  is 
reduced  to  a  minor  effect  by  choosing 


or,  equivalently, 


which  demands  a  very  large  number  of  sensors  or  telemetry  channels. 
This  may  not  be  convenient,  or  even  possible,  in  practice.  Such 
intervals,  as  remarked  earlier,  have  a  typical  length  scale  in  inches. 

The  influence  of  spatial  aliasing  is  illustrated  in  Fig.  4.1, 
where  the  sampled  dutterwortii  noise  spectrum  from  Fig.  2.3  is  shown 
for  several  uniform  sampling  intervals,  with  tne  postulated  signal  at 
vg.  The  wavenumber  spectrum  of  the  sampled  process  is  clearly 
periodic  with  regular  sensor  spacing.  As  a  tool  for  demonstrating 
effects  of  spatial  aliasing  noise,  fractional  bandwidth,  in  Fig.  4.1, 
has  been  purposely  exaggerated  compared  with  the  typical  spectrum  in 
Fig.  1.2.  Thus,  Fig.  1.2  gives  an  incomplete  picture  of  the  effective 
flow  noise  spectra  one  would  expect  with  discrete  arrays. 
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Sampled  (aliased)  spectra  in  Figs.  4.1b-c  are  obviously 
distortions  of  the  original  noise  spectrum  in  Fig.  4.1a.  Uniform 
sampling  causes  high  wavenumber  components  to  fold  down  onto  lower 
wavenumbers.  This  wavenumber  folding,  generated  strictly  by  spatial 
sampling,  clearly  increases  the  noise  level  and,  hence,  can  severely 
degrade  signal-to-noise  ratio  at  v  . 

Fig.  4. Id  shows  an  example  where  the  spatial  interval  is 
sufficiently  large  so  that  the  noise  is  nearly  uncorrelated  between 
sensors.  For  this  sampling  interval,  the  sum  in  v4.1)  scaled  by  v 
becomes  a  Rieman  sum,  which  represents  a  reasonable  approximation  to 
the  sampled  process  spectrum  level:  20ir(I  Is  ).  The  resultant 

K  K 

spectrum  for  the  sampled  process  is  nearly  wavenumber  white,  with  a 

spectrum  level  that  obviously  exceeds  the  peak  level,  21  /e  ,  for 

K  K 

continuous  observations  in  Fig.  4.1a.  Of  course,  the  power  spectrum, 
at  a  sensor,  obtained  by  integrating  (4.1)  over  any  periodic  interval, 
[v  ],  equals  I K ,  which,  from  our  remarks  in  Chapter  1,  increases 
witn  increased  tow  vessel  speed.  Conventional  beamforming  is  optimal 
at  this  spacing,  but  signal-to-noise  ratio,  at  the  beam  output,  is 
very  poor  compared  to  the  signal-to-noise  ratio  with  continuous 
observations.  A  common  practice  is  to  space  sensors  at  a  half 
acoustic  wavelength,  which  is  often  much  larger  than  2tt/&  .  Hence, 
detection  capability,  with  such  spacinys,  is  seriously  degraded 
compared  with  detection  capability  with  tighter  sensor  spacings. 
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One  infers  from  this  that  processing  with  discrete  observations 
always  results  in  unavoidable  spatial  aliasing,  making  it  impossible 
to  achieve  the  signal-to-noise  ratio  attainable  with  continuous 
observations. 

If  a  wide  band  of  frequencies  is  processed,  there  will  probably 
be  frequencies  at  which  the  noise  at  a  conventional  beamformer  output 
is  the  coherent  sum  of  noises  received  at  various  sensors.  High 
beamformer  output  noise  levels,  at  these  frequencies,  drastically 
degrade  detection  performance  unless  the  appropriate  postbeamformer 
filter  in  [2.89)  is  used. 

The  narrowband  version  of  the  sampling  theorem  [21,  pp.  34—3 b J , 
requiring  a  sampling  rate  of  2 tt / e ^ ,  is  irrelevant  here  because 
reconstruction  of  the  unsampled  noise  spectrum  cannot  be  accomplished 
from  raw  spatial  samples  taken  at  that  rate.  Once  discrete  sensors 
are  introduced,  spatial  aliasing  loss  has  happened  and  is  irreversible. 

4.3  Detector  Performance 

General  expressions  for  the  optimal  and  conventional  array  gains 
of  discrete  arrays  are  given  by  [2.55)  and  (2.69),  respectively. 
Expressions  for  the  detection  indices  are  given  by  (2.74),  (2.77),  and 
(2.79). 
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With  regular  sensor  spacing,  the  normalized  parameter, 
pk  =  pk^A^k’  (2.34),  completely  describes  the  noise  covariance 
matrix  Q^.  For  a  linear  array  of  M  regularly  spaced  sensors 
separated  by  a  distance  a,  the  exponential  form  of  (2.34)  leads  to  the 
Toeplitz  matrix: 
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(4.6) 


The  inverse  of  (4.6)  is  therefore  the  banded  matrix: 
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(4.7) 


which,  incidentally,  specifies  the  optimal  matrix  filter  shown  in 
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Fig.  4.2,  as  computed  from  (2.50).  The  filter  is  the  discrete  analog 
of  the  optimal  structure  for  continuous  observations  in  Fig.  3.1. 


The  optimal  and  conventional  array  gains  are  computed  (see 
Appendix  B)  in  a  straightforward  manner  from  (2.65)  and  (2.69): 
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\  =  e 


V  .  =  V,  -  V 

d  k  s 


Since  1/e,  is  the  noise  correlation  distance,  it  follows  that 


once  the  total  array  length  L  =  (M  -  1)a  exceeds  the  noise  correlation 
distance.  The  approximate  representation  of  conventional  array  gain 
is,  therefore. 
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(4.10) 


When  M  >«>,  (4.8)  and  (4.10)  approacn  the  same  limit  for  fixed  a: 

Gkc  «  Gko  “  -  _”| r|7  (‘  *  \K  I2  -  ?RA})  •  <4-u) 

Note  that  (4.11)  is  a  strong  function  of  u  as  long  as 
B^a  <<  1,  i.e.,  as  long  as  the  sensor  spacing  is  much  smaller  than 
the  noise  correlation  distance. 

More  interesting  than  the  limit  for  fixed  a  is  the  behavior  as  M 
varies  but  (M  -  1)a  =  L,  a  fixed  length  chosen  so  tnat  &.  L  >>  1.  To 

K 

deal  with  this  problem,  we  form  the  ratio  of  (4.8)  and  (4.10): 


(4.12) 
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After  some  algebraic  simpl ication,  this  reduces  to 


Re<V  - 


lA 


A  \2 


(1  -v 


(4-13) 


With  the  definition  of  from  (4.9),  (4.13)  appears  to  be  a 
function  of  the  four  parameters  8,  ,  v.  ,  a,  and  M.  However,  with 
the  constraints 

M  =  -  +  1  (4.14) 

A 


and 


=  constant, 

there  remain  only  two  independent  variables,  and  M.  for  any 
fixed  M,  the  ratio  GkQ/ Gkc,  considered  as  a  function  of  \>d, 
exhibits  a  maximum  L Gj^q/ Gkc max.  fig.  4.3a  plots  tnese  maxima  as 
a  function  of  M  for  various  values  of  8  L  >  3  (f>  ‘  <  U.02).  Once 

K  K 

8kL  exceeds  1  by  a  significant  factor,  the  ratio  of  array  gains 
barely  deviates  from  unity  for  all  possible  M,  hence  all  possible 
sensor  spacings. 
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It  is  relatively  easy  to  show  that  the  exact  expression  for 
optimal  array  gain  in  (4.8)  approaches  the  array  gain  with  continuous 
observations  in  (3.6)  in  the  limit  as  a  =  L/M  -  1  >  0  with  L  fixed. 


Inserting  the  definition  of  p,  in  (4.8),  we  obtain  the 

K 

following  equivalent  expression: 
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for  a  small  enough  so  tnat  a,  a  <<  1  and  av.  <<  1,  one  can 

K  0 

expand  the  exponentials  and  the  cosines  into  Taylor  series: 
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Then  performing  the  required  multiplication  and  collecting  terms,  we 
arrive  at 


iV'6kA  vd"  2 

G.  =  1  +  — =-  +  — t; - o  +  . ..  +  terms  in  a  and  higher  - 

k°  2  2  * 


Avd  3BkA 
-2Sk"  2 


(4.1bb) 


Taking  the  formal  limit  of  (4.15b)  as  a  =  L/M  >  0  with  L  fixed  and 
L  =  (M  -  1)a  «  iv1a,  we  ootain 


lim  Gko  " 
L/M  >  0 


BkL 
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kl 


w  +  1  • 


+  1 


(4.15c) 


which  is  identical  to  tne  result  in  (3.6). 

In  section  3.1,  we  observed  that  spatial  sampling  increases  the 
level  of  the  noise  wavenumber  spectrum  at  the  signal  wavenumber, 
v  .  The  increase  was  described  by  the  factor  (Fka(vs))/(Fk(vs)) 
in  (4.3).  We  now  demonstrate  oy  example  that  tne  same  factor 
cnaracteri zes  the  decrease  in  array  gain,  at  least  for  ekL  >>  1, 
i  .e. , 
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cont 


F<a(  V 


14.16) 


Fig.  4.3b  snows  two  curves:  (i)  uko  from  (4.3)  and  ^2)  | con t 

from  (3.6)  multiplied  by  the  ratio  (F, (vl)/(F,  (v  ))  from  (4.3). 

Both  are  plotted,  using  the  Corcos  model  for  e.  and  v,  in  (2.32) 
and  (2.33),  for  M  =  100,  a  =  0.01,  and  =  0.  The  two  curves  are 
indistinguishable  over  the  entire  range  shown,  indicating  that  the 
aliasing  noise  fully  accounts  for  tne  array  gain  degradation  as  one 
makes  the  transition  from  continuous  to  discrete  observations . 

4.3.1  Numerical  Examples 

Numerical  examples  can  now  be  computed  from  (4.8)  and  (4.9). 

Again  we  use  the  Corcos  model  values  for  3  and  v.  from  1 2 . 32 )  and 

K  K 

(2.33)  and,  in  the  interest  of  algebraic  simplicity,  we  assume  a 
broadside  signal  arrival.  No  loss  of  generality  is  implied  because 
one  can  always  choose  tne  origin  of  tne  wavenumber  space  at  the  signal 
wavenumber. 

Fig.  4.4a  shows  (4.11)  normalized  by  M  plotted  as  a  function  of 
to,  a / U  for  a  =  0.01,  i.e.,  a  fractional  bandwidth  of  0.01.  The 

K  C 
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low  array  gain  at  (w.  a/U  )  =  2nn,  n  =  1,  2,  is  due  to  the 

fact  that  at  these  points  the  noise  components  at  sensors  within  one 
correlation  distance  (1/s  )  of  each  other  are  very  nearly  in  phase. 
Hence,  subarrays  of  this  length  have  an  array  gain  close  to  one. 

Since  the  entire  array  can  be  decomposed  into  e^L  such  subarrays, 
the  total  array  gain  at  w^a/Uc  =  2nn  is  only  of  the  order  Bol¬ 
and  the  normalized  array  gain  of  Fig.  4.4a  is  GkQ/ M  «  B^a.  (More 
exact  calculations  from  (4.11)  yield  G^q / M  «  b^a/2)  .  It  is 
important  to  remember  that  poor  performance  at  wkA/Uc  =  2nn, 
n  =  1,  2,  ...,  results  from  the  fact  that  tne  noise  components  add 
coherently  at  these  frequencies  and  can  therefore  easi ly  dominate  the 
beamformer  output.  By  eliminating  tnese  undesirable  frequencies,  the 
modified  Eckart  filter  of  (2.89),  i.e.,  the  filter  optimized  for  i ts 
actual  input  signal  plus  noise  spectra,  achieves  dramatically  improved 
detector  performance. 


To  get  an  idea  of  the  extent  of  the  approximation  in  (4.11), 
normalized  but  exact  curves  from  (4.8)  and  (4.9)  are  plotted  in  Fig. 
4.4b  for  M  =  100  and  a  =  0.01.  Observe  that  once  WjA/U  >  it,  the 
difference  between  GkQ  and  G^c  quickly  becomes  trivial.  With 
a  =  0.01  and  M  =  100,  bkL  =  199/100) (a>k&/Uc)  •  Tdus,  the  tw0 
curves  coincide  once  b.,L  2l  More  detailed  examination  of  (4.8) 

k 

and  (4.9)  in  the  expanded  interval  [0-tt]  in  Fig.  4.4c  shows  that  Gkc 
actually  oscillates  about  the  curve  in  Fig.  4.4b  with  spatial 

frequency  Mv  and  thus  can  be  as  much  as  3-4  decibels  below  the 

k 
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average  in  the  vicinity  of  Au)k/Uc  =  u/10.  Conventional 
processing,  in  this  expanded  region,  is  obviously  closer  to  optimal 
when  Av^  is  an  odd  multiple  of  n/M.  Oscillations  decay  essentially 
to  zero  at 

u>kA/Uc  =  TT  /  2  . 

Because  of  the  approximate  equality  of  G,  and  G,  ,  the 

ko  KC 

detection  index  of  the  conventional  detector  using  the  modified  Eckart 

filter,  [see  (2.89)]  is  approximately  equal  to  dQ^ 

under  the  conditions  leading  to  (4.11).  On  the  other  hand,  if  one 
only  uses  the  simple  Eckart  filter,  v^/n  i.e>>  the  filter 

matched  to  the  i nput  signal  and  noise  spectra,  in  conjunction  with  a 
conventional  beamformer,  one  can  expect  drastically  degraded 
performance.  Specifically,  let 


y  (a  constant),  for  0  <_  a>k  £  w 
0  ,  elsewhere  . 


(4.17) 


Then,  from  (2.74)  and  (2.77), 
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where  N  =  Tw/2u  is  the  number  of  frequency  components  in  the  signal 
band. 

As  we  have  already  seen,  at  very  low  frequencies  the  3^L  >>  1 
assumption  would  be  violated  so  that  Gec  <  GeQ.  Hence,  the  ratio 
of  (4.18)  is  a  somewhat  conservative  estimate  of  the  performance  loss. 

Fig.  4.5  (uppermost  plot)  shows  (4.13)  plotted  as  a  function  of 
a,  using  array  gain  functions  approximated  by  4.11  with 

^  =  2w  .  (4.19) 

u 


Once  wa/Uc  _>  2ir,  the  ratio  of  detection  indices  becomes  insensitive 
to  the  precise  value  of  w,  as  one  might  expect  from  the  near-periodic 
structure  of  Fig.  4.4a.  Only  when  w  increases  sufficiently  so  that 
(a/Uc)wA  »  1  does  the  ratio  (4.18)  diminish  significantly  (because 
tne  high  frequency  noise  components  at  adjacent  sensors  are  now 
uncorrelated).  Thus,  array  gain  becomes  frequency  invariant  over  a 
major  portion  of  the  band. 
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More  exact  ratios,  also  shown  in  Fig.  4.5,  using  array  gain 
functions  from  (4.8)  and  (4.9),  show  clearly  that  the  approximation  in 
(4.18)  is  very  optimistic  for  noise  fractional  oandwidths  less  than 
0.01.  We  also  see  that  tne  detection  index  for  conventional 
oeamforming  with  modified  EcKart  filtering  is,  at  most,  b  decibels 
below  the  detection  index  with  optimal  array  processing.  The 
inaccuracy  of  (4.18)  becomes  unimportant  in  Fig.  4.5  when  a  increases 
to  0.1,  where,  presumably,  most  of  the  frequency  components  in 
Dandwidth  w  satisfy  the  large  condition. 

For  small  arrays,  the  approximations  leading  to  (4.11)  are 

inappropriate  and  one  must  again  deal  with  the  exact  (4.8)  and  (4.9). 

This  is  done  in  the  remaining  figures.  Array  gains  are  plotted  as 

functions  of  u>.  a/U  .  Wideband  detection  indices  are  based  on 
k  c 

spectra  of  the  form  defined  by  (4.17). 

Fig.  4.6a  shows  the  optimal  and  conventional  array  gains  for  a  2 
sensor  array  with  a  =  0.01.  At  a/J  =  (2n  +  l)u,  n  =  0,  1,  2, 

K  C 

...,  the  sensors  are  spaced  a  half  wavelength  apart  so  that 
conventional  beamforming  acnieves  optimal  noise  cancellation.  For 
full  wavelength  spacing,  the  noise  components  at  the  sensors  are  in 
phase,  just  as  the  signal  components,  and  no  signal  processor  can 
enhance  signal-to-noise  ratio.  For  these  reasons,  the  curves  coincide 
at  all  integral  multiples  of  n,  and  tne  array  gain  at  even  multiples 
of  it  is  0  decibels . 
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Fig.  4.6b  shows  the  wideband  detection  indices  tor  the  same  2 
sensor  array  and  fractional  bandwidth,  a  =  0.01,  plotted  as  a  function 
of  the  normalized  signal  bandwidth  w&/Uc.  The  upper  curve 
corresponds  to  the  optimum  detector,  the  middle  curve  to  tne  detector 
with  conventional  beamforming  followed  by  modified  Eckart  filter,  and 
the  lower  curve  to  the  detector  with  conventional  beamforming  followed 
by  the  simple  Ec<art  filter.  For 


tne  modified  Eckart  filter  yields  very  little  performance  improvement 
and  the  loss  relative  to  optimal  performance  is  serious  (22-30  decibels) 
However,  once 


IT  >  11  ’ 

c 

the  conventional  detector  with  a  modified  Eckart  filter  comes  within 
9  decibels  of  the  optimal  detector,  whereas  the  conventional  detector 
with  a  simple  Eckart  filter  falls  25-31  decibels  below  the  optimum. 

In  this  range,  the  separation  of  the  three  curves  does  not  vary 
drastically,  indicating  tnat  relative  performance  is  insensitive  to 
the  specific  values  of  w  and  Uc. 
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Figs.  4.7-4.10  snow  the  dependence  of  array  gain  on  a  for  a  fixed 
2  sensor  array.  With  a  =  1,  the  noise  is  near-white  and  both 
conventional  and  optimal  processors  have  an  array  gain  close  to 
2  decibels.  As  the  fractional  bandwidth  of  the  noise  decreases,  the 
maximum  array  gain  increases  sharply.  The  minimum  must,  of  course, 
remain  0  decibels  for  the  reasons  mentioned  in  the  discussion  of 
Fig.  4.6a. 

Fig.  4.11  shows  the  dependence  of  the  wideband  indices  on  a.  The 

array  gains  required  for  the  computation  were  taken  from  Figs. 

4.7-4.10  and  the  normalized  signal  bandwidth  was  held  fixed  at 

wa/U  =  2 it .  Once  a  is  small  enough  so  tnat  the  difference  between 
c  3 

2  2 

d  and  d  i  c  ,  ,  becomes  significant,  use  of  the  modified  Eckart 

0  C  I  tCK dr  t 

filter  recovers  a  major  fraction  of  the  loss. 

The  final  series  of  curves  deals  with  the  effect  of  the  number  of 
sensors,  M,  on  relative  performance.  Fig.  4.12  snows  the  array  gains 
for  3-  and  5-element  arrays  with  a  =  0.01.  Conventional  and  optimal 
array  gains  are  approximately  equal  near  the  M  -  1  frequencies 
satisfying 

^  ^  p  ,  p  =  1,  2,  ...  (M  -  1)  .  (4.20) 

c 

Fig.  4.13  gives  trie  wideband  detection  indices  normalized  with 
respect  to  tne  assumed  constant  signal-to-noise  ratio  in  (4.17)  and  M 
as  functions  of  M  for  wa/Uc  =  2n.  Once  again,  use  of  the  modified 
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EcKart  filter  recovers  much  of  the  loss  incurred  in  using  conventional 
beamforming  with  a  simple  EcKart  filter.  The  extent  of  that  recovery 
grows  with  M,  becoming  complete  when  BK(M  -  1)a  =  BKL  >>  1  so  that 
performance  is  characterized  by  Fig.  (4.5). 

4.4  Summary 


We  examined  spatial  aliasing  with  arrays  of  uniformly  spaced 
sensors  and  found  aliasing  to  De  the  major  cause  of  detection  loss 
attendant  on  discrete  observations.  Attendant  loss  in  array  gain 
relative  to  array  gain  with  continuous  observations  is  directly 
traceable  to  the  increased  noise  level,  due  to  aliasing,  at  the  signal 
wavenumber. 

When  s,L  »  1,  a  condition  usually  satisfied  by  towed  arrays  in 
flow  noise,  conventional  array  gain  is  essentially  optimal  with 
uniform  sensor  spacing,  detection  performance,  in  this  case,  is, 
however,  limited  by  spatial  aliasing,  which  must  degrade  the 
signal-to-noise  ratio.  Thus,  signal  detectaoi 1 ity  is  degraded 
relative  to  the  case  with  space-continuous  data.  Spatial  aliasing  is 
avoided  only  if  spatial  observations  are  effectively  continuous,  with 
extremely  small  sensor  intervals  so  that  discrete  observations  pass  in 
the  limit  to  continuous  oDservations.  For  the  flow  noise  problem, 
sensor  spacings  on  the  order  of  inches  are  required  to  approach  the 
continuous  mode  of  operations. 
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As  with  continuous  observations  in  Chapter  3,  use  of  a  relatively 
simple  modification  of  the  Ecxart  filter  can  very  substantially 
enhance  the  performance  of  a  detector  using  conventional  beamforming, 
making  it  more  nearly  competitive  with  the  optimal  detector.  Major 
gains  can,  of  course,  be  made  only  when  a  reasonably  wide  band  of 
frequencies  is  processed  and  tne  noise  has  sufficient  spatial 
structure  so  that  optimal  and  conventional  array  gains  vary 
significantly  with  frequency.  The  flow  noise  associated  with  towed 
arrays  has  the  general  form  of  the  model  used  in  the  examples.  While 
its  relative  bandwidth  need  not  oe  as  low  as  0.01,  it  would  probably 
be  sufficiently  concentrated  in  wavenumber  space  to  make  the  modified 
Eckart  filter  attractive  for  wideband  detection. 

A  clear  drawback  of  the  modified  Eckart  filter  is  its  dependence 
on  the  (conventional)  array  gain  and  hence  on  the  spatial  structure  of 
the  noise  field  (except  for  the  continuous  line  array  discussed  in 
Chapter  3).  This  structure  is  not  likely  to  be  known  in  advance,  so 
tnat  some  form  of  adaptation  is  required. 

Here  again,  one  needs  only  a  single  adaptive  filter  rather  than 
the  set  of  M  filters  required  for  optimal  processing  of  data  from  M 
sensors.  For  flow  noise  encountered  by  towed  arrays,  the  primary 
unknown  parameter  would  probably  be  Uc,  a  quantity  closely  related 
to  tow  speed.  Adaptation  could,  therefore,  be  a  relatively  simple 
matter;  one  could,  for  example,  estimate  the  noise  spectrum  at  the 
conventional  beamformer  output  by  simply  estimating  U  . 
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Fig.  4.1 .  Wavenumber  Spectral  Overlap 
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SUMMERS  MULTIPLIERS 


Fig.  4.2.  Optimal  Linear  Combiner  for  M  Equally  Spaced  Sensors 
in  [-L/2,  L/2]  With  Butterworth  Noise  Model 
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Fig.  4.3a  Maximum  Ratio  of  Optimal  to  Conventional 
Array  Gain 
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Fig.  4.3b  Optimal  Array  Gain  With  Spatial 
Aliasing  Corrections  on  Optimal  Array  Gain  With 
Continuous  Observations 
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Fig.  4.4.  Optimal  and  Conventional  Array  Gain  Versus 
tokA  U.  for  a  =  0.01.  L  =  (M-l)A 
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Fig.  4.5.  Wideband  Detection  Index  Advantage  With  Modified 
Eckart  Filtering  Over  Simple  Eckart  Filtering 
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Fig  4  6a  Optimal  and  Conventional  Array  Gain  for  a  2  Sensor  Array 

With  a  =  0.01 


Fig.  4.6b.  Wideband  Detection  Indices  for  2  Sensors  With  a  =  0.01 
Plotted  as  a  Function  of  Normalized  Bandwidth,  wA 

Uc 
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Fig.  4.7. 
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Fig.  4.7.  -  4.10.  Array  Gain  Dependence  on  a  =  falw  for  a 
2  Sensor  Array,  M  =  2 
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Fig.  4.12.  Array  Gain  for  3  and  5  Sensors,  a  =  /3i  n  =  0.01 
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NUMBER  OF  SENSORS  M 

Fig.  4.13.  Normalized  Detection  Indices  Versus  Number  of  Sensors, 
M  With  a  =  i\  =  0-01  and  wA,  U,  =  2 n 
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CHAPTER  FIVE 


ARRAY  DESIGN  CONSIDERATIONS 


5.1  Introduction 

In  Chapters  3  and  4  array  geometry  was  given.  We  must  now  deal 
with  criteria  that  affect  the  choice  of  geometry  in  practice.  In  the 
continuous  case,  the  primary  constraint  is  that  the  introduction  of 
continuous  phase  shifts  is  impractical  so  tnat  one  must  ideally  work 
with  sensor  lengtns  small  compared  with  acoustic  wavelengths. 

In  the  discrete  case,  the  primary  constraint  is  a  limitation  on 
the  number  (M)  of  available  sensors.  The  number  must  often  be  small 
enough  so  that  uniform  spacing  over  the  available  aperture  results  in 
the  flow  noise  being  independent  from  sensor  to  sensor,  yielding  an 
array  gain  of  only  M.  By  using  groups  of  tightly  spaced  sensors,  one 
can  get  a  much  nigher  array  gain  at  certain  frequencies  and  hence 
achieve  greatly  improved  detector  performance  (either  optimal  or 
conventional ) . 

5.2  Space-Continuous  Sensors 

The  first  detector  structure  we  consider  is  shown  in  Fig.  5.1. 

It  consists  of  conventional  beamformi ng  on  N^  space-continuous 
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uniformly  sensitive  sensors  (extended  sensors)  of  length  1.  Spatial 
aliasing  is  thus  avoided  at  the  extended  sensor  level  out,  of  course, 
not  at  the  beamformer  input  level.  Extended  sensors  are  separated  by 
distance  D,  which  we  assume  to  be  sufficient  for  the  noise  to  be 
statistically  independent  oetween  extended  sensor  outputs.  Hence, 
conventional  beamforming,  which  simply  steers  extended  sensor  outputs 
to  align  signal  components,  is  formally  optimal. 

For  most  of  the  discussion,  equally  spaced  subapertures  are 
assumed.  Tnis  assumption  is  unnecessary.  It  is  simply  convenient. 

AM  of  the  arguments  to  follow  remain  valid  for  irregularly  spaced 
subapertures .  The  main  requirement  is  that  tne  outputs  of  subapertures 
are  uncorrelated  under  tne  noise-only  hypothesis.  It  is  important  to 
keep  in  mind  that,  if  the  noise  is  uncorrelated  oetween  subarrays,  the 
subarrays  may  oe  distributed  in  tne  observation  interval  in  any  manner 
whatever,  leaving  the  results  unchanged. 

In  passing,  we  emphasize  tnat  one  is  driven  to  this  type  of  array 
construction  because  of  constraints  on  the  number  of  sensors. 

Otherwise,  one  could,  in  principle,  fully  populate  the  observation 
interval  with  a  large  number  of  equally  spaced  extended  sensors  of 
very  small  length  and  thus  achieve  near  optimal  performance  with 
simple  conventional  beamforming  as  in  Chapters  3  and  4,  as  long  as  the 
overall  product  is  sufficiently  larger  tnan  unity. 
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Appendix  B  derives  the  overall  array  gain  for  tne  structure  shown 
in  Fig.  5.1  as  a  direct  extension  of  the  results  obtained  in  Chapter  3. 
The  array  length  is  now  1,  but  no  steering  is  possible  within  the 
space-continuous  extended  sensor.  Hence, 

G„r  -  N,  $  2  G.  i  ,  (5.1) 

Kr  1  c  kc|suD  c 

where 


G 


kc 


sub  c 


A 


1 


COS  v^l  + 


2Vk 


A  ‘  \ 


2  \ 


(5.2) 


is  the  array  gain  for  the  extended  sensor.  There  snould  be  no 
confusion  between  (5.1)  and  (5.2).  The  quantity,  Sc  G^c  .  in  (5.1) 
is  the  unsteered  array  gain.  Equation  (5.2)  gives  the  array  for  a 
broadside-looking  direction.  Excepc  for  the  replacement  of  by 
(no  steering),  (5.2)  is  analogous  to  the  array  gain  for  the 
entire  (steered)  array  in  Chapter  3  (3.14)  and  is  common  to  all 
channels.  Therefore, 


S 


c 


2 


(5.3) 
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is  the  wavenumber  response  of  the  extended  sensor  to  signal.  is 
the  array  gain  factor  from  conventional  beamforming  on  the 
extended  sensors,  witn  noise  between  sensor  outputs  statistically 
independent. 

The  principles  and  reasons  behind  using  modified  Eckart 
filtering,  rather  than  simple  Ecxart  filtering,  in  Fig.  5.1  are 
exactly  the  same  as  in  Chapter  3.  The  modified  Eckart 
(postbeamformer)  filter,  C^,  from  (2.89)  should  be  chosen  according 
to  its  input  signal  and  noise  spectra,  not  the  signal  and  noise 

spectra  at  each  sensor. 

Let  us  now  consider  that  the  signal-to-noise  ratio  is  constant 
over  the  signal  band: 


S/I  (a  constant)  ,  k  =  1,  2,  ...»  N 
o,  elsewhere  , 

N  =  Tw/2tt  . 


(5.4) 


Then,  using  conventional  beamforming  oil  the  extended  sensor  outputs 
(H  =  V  in  Fig.  5.1)  and  simple  Eckart  filtering,  =v/S^/Ij<, 

— k  "rs. 

we  obtain  from  (2.77)  and  Appendix  6 
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(5.5) 


From  (5.3)  and  the  properties  of  the  "sine"  function, 


Hence, 


(5.5b) 


Equality  applies  in  (o.5a)  and  (5.5b)  only  when  o$  =  ir/2.  In  such  a 
case,  the  signal  arrives  broadside  to  the  array  axis.  Hence,  the 
conventional  detection  index  in  (5.5)  is  degraded  relative  to  a  system 
using  conventional  beamforining  with  steering  within  the  extended 
sensors.  Tne  degradation  is  clearly  large  if  t tie  signal  bandwidtn  is 
concentrated  largely  outside  tne  mainlobe  response  in  (5.3). 
Conversely,  one  nearly  avoids  this  aegradation  if  the  largest  signal 
frequency  is  well  inside  the  mainlobe  of  (5.3).  For  tnis  latter  case, 
tiie  extended  sensor  length  is  muen  smaller  than  any  signal  wavelength 
along  the  array. 


TR  7997 


Using  a  common  modified  Eckart  filter,  -V^S^/I,^  G^y, 


in  conjunction  with  the  conventional  beamformer,  in  Fig.  5.1,  we 
obtain  from  (2.79),  (5.1),  and  (5.4) 


c  modified 
Eckart 


(5.6) 


As  in  (5.5),  performance  in  (5.6)  is  degraded  relative  to  performance 
with  steering  within  extended  sensors,  insofar  as  the  signal  bandwidth 
is  outside  the  mainlobe  of  (5.3).  In  many  ways,  the  discussion  is  the 
same  here  as  in  Cnapter  3.  The  idea  of  not  steering  within  extended 
sensors,  which  causes  signal  power  suppression  prior  to  conventional 
beamforming,  is  tne  main  difference. 

Here,  as  in  Chapter  3,  we  might  expect  conventional  beamforming 
followed  by  modified  Eckart  filtering  to  recover  a  substantial  part  of 
the  performance  loss  incurred  with  simple  Eckart  filtering  and 
conventional  beamforming.  This  turns  out  to  oe  false  in  some  cases, 
e.g.,  when  Sc2  «  1  throughout  the  processed  frequency  band.  Note  that 
the  best  frequency  filter,  in  this  case,  requires  the  weight  of  the 
common  subaperture  gain,  including  in  (5.3),  in  addition  to  the 
simple  Eckart  f i Iter,  VSK/ 1^. 
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From  (5.1)  we  see  tnat  G^y  varies  directly  with  the  extended 
sensor  broadside  array  gain  in  (5.2)  and  extended  sensor  response  to 
signal.  Therefore,  the  sum  in  tne  denominator  of  (5.5)  is  affected 
most  heavily  by  frequencies  wnere  extended  sensor  array  gain  is  small, 
e.g.,  at  frequencies  satisfying 


=  (2n  +  1 ) 7i  ,  n  =  0,  1,  2,  .... 


(5.7) 


At  these  frequencies,  1  is  an  odd  multiple  of  tne  noise  half 
wavelength,  evaluated  at  u>k« 


Tne  sum  in  (5.6)  is  influenced  most  heavily  by  frequencies  where 


G)<c  |  is  large  and  relatively  little  by  frequencies  satisfying 


|SUQ  C 

(5.7).  Simple  Eckart  filtering  is,  once  again,  poorly  matched  to 


conventional  beamforming  if  there  is  significant  extended  sensor  array 


? 

gain  variation  with  frequency,  quite  apart  from  the  effect  of  Sc  in 
(5.3).  Arguments  made  in  Chapter  3  in  connection  with  steered 
aperture  thus  provide  us  with  a  basis  for  discussion  here. 


Consider  first  that  extended  sensor  length  is  much  smaller  than  any 
signal  wavelength,  i.e.,  v^/2  «  1  in  (5.3)  for  k  =  1,  2,  3,  ...,  N. 
Then  from  (5.5)  and  (5.6), 
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d  I 

c  [modified  EcKart  _  l 

d  2 1  N 

c  |simple  Eckart 

This  is  the  same  for  tne  entire  space-continuous  array  as  for  a  simple 
extended  sensor.  We,  therefore,  simply  appeal  directly  to  Chapter  3 
for  all  relevant  conclusions. 

Conventional  processing  with  modified  Eckart  filtering  is  better 
than  conventional  processing  with  simple  Eckart  filtering  unless  tne 
extended  sensor  array  gain  is  frequency  independent.  The  extent  of 
the  improvement,  previously  discussed  in  connection  with  tne  numerical 
examples  of  Fig.  3.4,  is  suostantial,  especially  for  small  akl 
products  and  small  noise  fractional  oandwidths, 

Hence,  if  one  constructs  an  array  of  extended  sensors  that  are 
simultaneously  small  compared  with  signal  wavelengtns  and  with  the 
noise  spatial  correlation  length,  then  modified  Ecxart  filtering  can 
potentially  enhance  performance  of  the  conventional  detector,  making 
it  more  competitive  with  the  optimal  detector.  Here  again,  as  in 
Chapter  3,  the  drawDack  of  modified  Eckart  filtering  is  the  strong 
dependence  on  the  extended  sensor  array  gain  in  (5.2)  and  hence  on  the 
spatial  structure  of  the  noise  field,  which  usually  varies  with 
conditions.  In  such  cases,  one  generally  needs  to  estimate  the 
spectral  parameters  vk  and  eR  in  order  to  construct  the  modified 


E  E 


N  P  2 
bkc 


sub  c 


:=1  j-1  Sjc 


2 


(5.8) 


sub  c 
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Eckart  filter  from  (5.2):  remember,  for  the  Corcos  flow  noise  model, 
convective  noise  speed,  l)  ,  is  sufficient  for  this  construction. 

We  next  consider  cases  for  detection  performance  witn 

? 

conventional  Deamforming  when  Sc  in  (5.3)  is  not  equal  to  1  over 
the  entire  signal  Dand.  We  shall  confine  this  discussion  to  the  large 
6.  1  regime,  where  extended  sensor  response  to  signal  tends  to  be 
more  of  a  problem. 

We  ask  if  there  is  substantial  detection  enhancement  possiole 
with  modified  Ecxart  filtering  at  the  oeamformer  output  when  there  is 
no  steering,  within  extended  sensors,  to  align  tne  signal  properly? 
When  the  answer  is  yes,  we  ask  if  tins  enhancement  is  about  the  same 
as  the  enhancement  accomplished,  under  the  same  conditions,  with 
steering  within  extended  sensors?  If  so,  tnen  why  bother  to  steer? 
Performance  loss  due  to  signal  suppression  at  tne  extended  sensor 
level  could  then  be  recovered  at  the  beam  output  with  appropriate 
scalar  filtering.  It  turns  out  that,  in  general,  the  answer  to  the 
initial  question  is  no. 

We  show  tnis  using  our  usual  noise  model  witn  g,  and  v,  given 
by  (2.32)  and  (2.33)  witn  spectral  functions  specified  at  eacn 
frequency  by  (3.29).  Let  us  also  assume  a  sufficiently  large 
T(u>2  -  w^)  product  so  that  sums  in  (5.5)  and  (5.6)  can  be 
approximated  by  integrals.  From  (5.5)  and  (5.6),  one  obtains 
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=  I_  N  2  Si  Jl_ 

modified  Eckart  ^  ^  1^  4Uc^ 


and 


j_  N  2  £  _Li 

simple  Eckart  ^  ^  1^  4U  ^ 


dw 


(li  o 

A 

*U)^  0) 


do> 


Bkl  »  1  ,  (5.10) 


which,  again,  are  simple  extensions  of  our  results  in  Chapter  3. 
Making  the  changes  of  variaoles 


U>  1 

u  -  2^  cos  es 


du  =  cos  &s  du, 


u,  = 


(*)^  I 

2c 


COS  9, 


(5.11) 


u„  = 


IE  cos  es 


with  c  =  speed  of  sound,  we  have  the  ratio 
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[modified  Eckart 


|simple  Eckart 


no  steering 


f2 

- 

V  u 


du 

2 


(5.12) 


Again,  when  3 k 1  >>  1  is  violated,  (5.12)  is  an  optimistic  estimate 
of  ennancement.  Note  that  (5.12)  is  the  same  as  (3.32)  when  the 
signal  arrives  Droadside  to  the  array,  i.e.,  e$  =  n/2  so  that 
2 

Sc  in  (5.3)  becomes  1  in  trie  signal  band. 

Certainly,  enhancement  without  steering  within  extended  sensors 
as  given  by  (5.12)  must  always  be  greater  tnan  1.  This  follows 
directly  from  (5.12)  by  using  tne  Schwartz  inequality: 


0<U^<U^<00. 


(5.13) 
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Therefore, 


(5.14) 


Hence,  (5.12)  exceeds  1  for  aroitrary  nonzero  intervals  [u^,  u^]. 
Equality  applies  in  (5.14)  when  »  u^. 


Under 


the  steered  condition,  all  S 1  -  -  terms  in  (5.12)  go  to 


unity  with  the  result 


{modified  Eckart 


simple  Eckart 


|  2  Au  f 


du 

2 


steered 


u^  u 


f 


du 


u  = 


2c  5 


V  »  1  • 


.4  . 


3  \  u 


(5.15) 


whi 


ich  is,  of  course,  identical  with  (3.32). 


We  now  want  to  show  that  failure  to  steer  can  sometimes  result  in 
serious  performance  loss.  To  demonstrate  this,  we  have  plotted  the 
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exact  equations  in  (5.12)  and  (5.15)  as  functions  of  u^  for  three 
fixed  values  of  u^.  Results  are  shown  in  Fig.  5.2. 

Tne  results  indicate  tne  following  main  points.  Performance  gain 
provided  by  the  modified  Ecxart  filter  in  conjunction  with 
conventional  beamforming  on  extended  sensors  is  nearly  the  same 
whether  or  not  one  steers  within  extended  sensors  if 


u 


2  ~ 


2 

2c  C0S  es 


(5.16) 


Observe  that  each  curve  describing  enhancement  witn  steering  departs 
dramatically  from  the  corresponding  curve  for  the  nonsteered  condition 
at  tt/2,  which,  incidentally,  is  where  the  signal  response  in  (5.3) 
drops  to  3  decibels  below  its  value  at  v^/2  =  0. 


If  the  condition  in  (5.16)  is  satisfied,  tne  extended  sensor 
length  is  small  compared  to  signal  wavelengths.  Fig.  5.2  indicates 
that  one  then  achieves  approximately  5.12  decibels  of  conventional 
detection  enhancement  when  u^  =  tt/40,  using  the  simple  post- 
beamformer  filter  in  (3.33)  consisting  of  a  voltage  gain  increasing 
linearly  with  frequency  in  the  frequency  interval  specified  by  (5.16). 

Fig.  5.2  also  suggests  that  relative  performance  is  best  when  the 
normalized  lower  frequency  is  smallest: 
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— 7T-  cos  9  «  1  .  (5.17) 

2c  s 

Conditions  (5.16)  and  (5.17),  togetner ,  suggest  the  following 
procedure  for  determining  suitable  lengtns  for  extended  sensors. 

Extended  sensor  length  1  should  be  no  longer  than  one-fourth  of  the 
smallest  effective  signal  wavelength  =  xs/cos  es), 

evaluated  at 

Observe  tnat  there  is  insignificant  gain  provided  by  the  modified 
Ecxart  filter  when  u1  =  W4.  This  is,  of  course,  to  be 

expected  since  then  Sc^  is  small  throughout  most  of  the  processed 
signal  band.  One  can  easily  generalize  this  to  the  condition 

u,  >  1  . 

For  instance,  with  u^  =  2n  and  u2  =  4tt  inserted  into  (5.12),  one 
finds  that 

? 

d 

c  [modified  Ecxart  -1.52  (5#18) 

.  2 

d  . 

1  simp le  Eckar  t  |  n0  steer-jng 

This  confirms  that  the  modified  Eckart  filter  construction  in  (5.1)  is 
not  used  efficiently  in  this  case. 
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For  fixed 

up  to 


and  incident  angle  e$,  we  can  process  frequencies 


ttC 

“2  =  21  cos  es 


(5.19) 


and  still  De  essentially  compatible  with  conditions  needed  to  use  the 
relatively  simple  voltage  gain  postbeamformer  filter  in  (3.33). 
Clearly  one  gains  little  by  processing  frequencies  any  higher  than 
(5.19)  with  only  this  filter. 


The  main  effect  of  the  signal  incident  angle  is  to  increase  the 
upper  frequency,  as  ©s  approaches  n/2,  i.e.,  as  one  approaches 
broadside  signal  incidence.  It  is  notable  that,  even  for  endfire 
incidence  (e.g.,  e$  =  0),  the  quantity  c/1  is  sufficiently  large  in 
most  towed  array  applications  so  tnat  is  usually  of  the  order  of 
nundreds  of  hertz. 


The  other  region  of  interest  in  Fig.  5.2  is  where  the  normalized 
frequency  (w~l/2c)  cos  ©s  exceeds  W2,  the  point  where  the  signal 
response  in  (5.3)  is  down  by  3  decibels  from  its  value  at  v  , /2  =  0. 

For  the  case  u^  <  1,  curves  representing  performance  recovery 
for  the  unsteered  condition,  in  Fig.  5.2,  approach  nearly  constant 
values  as  u^  exceeds  u/2.  On  the  other  hand,  curves  representing 
recovery  for  tne  steered  condition  remain  monotone,  increasing  with 
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u  The  additional  noise  admitted  by  the  simple  Eckart  filter  makes 
the  difference.  With  simple  Eckart  filtering  following  the 
beamformer,  we  are  admitting  less  ano  less  signal  but  more  and  moi  e 
noise  as  Up  increases  beyond  n / 2 .  In  contrast,  the  band  limitation 

of  S  2  in  the  Eckart  filter,  admits  less  signal  and  noise, 
c 

It  is  easy  to  verify  tnat  tne  quotient  of  forms  containing 
sin  u/u  terms  in  (5.12)  are  less  than  unity  and  nearly  invariant  with 
Up  when  Up  exceeds  w / 2 .  Equation  (5.12)  then  reveals  that 
recovery  with  modified  Eckart  filtering  is  traceable  directly  to 
signal  frequencies  near  u^,  where  the  large  e^l  product  array 
gain  in  (5.2)  is  lowest  for  the  Droadside-looking  direction;  hence, 
the  noise  level  admitted  by  the  simple  Eckart  filter  is  highest.  To 
snow  this,  consider  the  case 

U1  =  lo  and  U2  >  2 

Then,  from  (5.12),  we  have 


«  0.35  ,  (5.20) 

Up  >  it/  2  . 


r 

it/ 40 


2-4  0 

sin  u  2  . 
- , —  u  du 


u 


1 


.  2 

sin  u 


du 


/  40 
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Hence, 


d 


d 


c 


c 


2 


2 


modified  Eckart 
simple  Eckart 


no  steering 


0.35 


f  1 

40  1_ 

|  u2  du  =  0.35 

IT  Up 

tt/40 


=  4.46  (6.5  decibels) , 


U2  >  W  2 , 


(5.21) 


which,  as  in  Fig.  5.2,  is  obviously  invariant  with  Up,  under  the 
specified  conditions. 

A  similar  argument  can  De  made  for  the  other  examples  in 
Fig.  5.2.  With  u^  =  n/10  and  Up  >  it/' 2,  the  ratio  in  (5.20) 
increases  to  approximately  0.49:  tne  detector  witn  simple  Eckart 
filtering  also  encounters  a  proportiona  1  ly  smaller  noise  level  at  the 
low  edge  of  the  processed  band  with  u^  increased  from  tt/40  to  tt/10. 
This,  then,  results  in  tne  relatively  small  performance  recovery  of 
only  1.98  decibels  compared  with  the  u^  =  tt/40  case. 

We  visualize  that  substantial  signal  energy  is  rejected  by  the 
extended  sensor  when  Up  exceeds  tt/2.  Processing  normalized 
frequencies  beyond  tt/2  is  probably  not  worthwhile  for  two  obvious 
reasons.  First,  one  requires  a  relatively  complicated  postbeamformer 


135 


TR  7997 


filter.  Secondly,  one  improves  performance  only  si i gut ly  with  such  a 
filter  compared  with  tne  improvement  made  when  processing  up  to  only 
ir/4.  We  cannot  over  emphasize  the  fact  that  performance,  in  the 
latter  interval,  is  achieved  with  the  modified  Eckart  filter  in 
(3.33),  consisting  of  a  relatively  simple  voltage  gain  increasing 
linearly  with  frequency  and  independent  of  specific  operating 
conditions.  That  filter,  coupled  with  a  low-pass  filter  cutting  off 
at  u2  =  tt/4 ,  should  therefore  represent  a  practically  reasonable 
compromise. 

There  are  certain  narrow  frequency  regions  where  the  extended 
sensor  signal  response  in  (5.3)  introduces  large  variations  in  beam 
output  signal  power,  e.g.,  in  the  vicinity  of  it  wnere  there  are 
relatively  large  changes  in  the  extended  sensor  response  with 
frequency.  Hence,  performance  with  no  steering  applied  within 
extended  sensors  is  generally  improved  over  performance  witn  steering 
within  sensors  in  such  regions,  but  only  by  a  small  amount. 

For  instance,  with  u^  =  it  and  u^  =  3tt / 4 ,  we  have  by  direct 
numerical  evaluation  of  (5.12)  and  (5.15)  tnat 

=  1.58  .  (5.22) 

no  steering 


u  modified  tcKart 


simple  Eckart 
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and 


modified  Eckart 


simple  Eckart 


steered 


1.03  . 

(5.23) 


We  conclude  from  the  above  that  the  modified  Eckart  filter  recovers 
some  of  the  losses  due  to  high  frequency  signal  rejection.  However, 

even  in  narrowband  frequency  regions  of  maximum  variation  of 

? 

S  ,  tne  recovery  is  so  small  as  to  make  construction  of  the  modified 
c  J 

Eckart  filter  in  (5.6)  unattractive. 


When  the  e^l  product  is  sufficiently  larger  than  1  and  the 
unused  aperture  space  in  Fig.  5.1  is  a  small  fraction  of  the  available 
observation  interval  L,  then 


N,1  =  L  . 

It  follows  from  (5.6)  with  v  1/2  <<  1  tnat 


d  ^  I  modified  =  /  Id, 
c  |  Eckart  Y' 

is  nearly  the  optimal  detection  index  for  interval  L  computed  in 
Chapter  3. 


■f  ■£ 


(5.24) 
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5.3  Discrete  Sensor  Subarrays 

The  constraints  here  are  very  different.  There  is  no  longer  any 
reason  to  demand 


1  «  x $  (5.25) 

and  optimal  processing  within  each  subarray  is  at  least  potentially 
possible.  Also, 


1  »  l/s„  (5-26) 

is  no  longer  a  natural  condition,  particularly  with  small  numbers  of 
sensors  in  the  subarrays.  Tne  major  constraint  is  that  the  total 
number  of  sensors  ( M ) ,  limited  by  practical  considerations,  is  too 
small  to  fill  the  available  aperture.  Clearly,  with  unlimited  sensors 
and  channels,  we  could  fully  populate  the  array  and  achieve  near 
optimal  performance  with  conventional  processing  on  individual  equally 
spaced  sensors  so  long  as  >>  1.  What  we  are  trying  to  do  here 
is  to  get  the  noise  cancellation  and  reduced  spatial  aliasing 
available  with  small  sensor  spacings  while  retaining  large  overall 
aperture  lengths  for  tracking  purposes. 
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Tne  general  overall  array  gain  for  this  detector  structure  with 
conventionally  processed  clusters  (subarrays) (see  Fig.  5.3),  derived 
in  Appendix  B  and  applicable  to  irregularly  spaced  subarrays,  is 


1  sub 


(5.27) 


i=l 


Equation  (5.27)  is  analogous  to  tne  array  gain  for  the  entire  array  in 
(4.9),  Chapter  4.  The  basic  discussion,  numerical  results,  and 
conclusions  there  are  tnus  relevant  nere  (see  Figs.  (4.6a),  4.7,  4.10, 
and  4.12).  The  detection  index  for  the  detector  structure  in  Fig.  5.3 
is  obtained  from  (2.86),  with  array  gain  from  (5.27): 


(5.28) 
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Except  for  tne  factor  ,  this  is  analogous  to  the  detection 
performance  with  full  array  processing  in  Chapter  4.  Hence,  the 
analogous  expression  for  the  relative  wideband  performance  gain  for 
the  conventional  detector  with  a  modified  Eckart  filter  over  the 
conventional  detector  with  simple  Eckart  filtering  serves  as  a  useful 
measure  of  improvement  here.  With  signal  and  noise  spectra  satisfying 
(5.4), 


c  [modified 
I Eckart 


[simple 
I Eckart 


(5.29) 


Performance  examples  with  tne  Corcos  model  shown  in  Fig.  4.5, 

4.6b,  4.11,  and  4.13  are  directly  applicable  to  (5.29).  The  modified 

Eckart  filter  recovers  much  of  the  loss  incurred  when  simple  Eckart 

filtering  follows  conventional  beamforming  on  subarrays  for  the  same 

reasons  as  in  Chapter  4.  Substantial  recovery  is  possible  when  a 

reasonable  wide  band  of  frequencies  is  processed  and  there  are  large 

variations  in  subarray  array  gain  with  frequency.  As  in  Fig.  4.13, 

where  detection  indices  are  normalized  by  the  square  of  the  product  of 

„  2 

the  total  number  sensors  used  times  S/I  in  (5.4),  i.e.,  (N^  S/IJ  , 
the  extent  of  the  recovery  increases  with  M^. 
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Analogous  to  processing  over  the  full  array,  one  can  improve 
detection  performance,  especially  in  the  small  6^,1  regime,  with  the 
detector  structure  in  Fig.  5.3,  by  using  optimal  processing  within 
clusters  (see  Fig.  4.2).  Such  an  instrumentation  involves  complex 
weighting  and  phase  shifting  of  sensor  outputs  and  results  in  an 
expression  for  overall  array  gain  given  by 


GkD  "  N1  Gko 


¥i 


sub 


1  -  e 


9  -2ao>.  A/U 

-2“V/Uc  11  +  11 


,  „  1_,  -«”k1/uc  cos  (Vj)  , 

L  1  M/e  *  U 

1  c 


\>S  <<  , 


(5.30) 


which  differs  from  the  expression  for  full  array  processing  in  (4.8) 
only  oy  the  appearance  of  in  place  of  M.  This  expression 

then  can  be  applied  directly  to  (2.79)  to  obtain  the  detection  index 
for  the  detector  structure,  as  in  Fig.  5.3,  with  conventional 
beamforming  on  optimally  processed  cluster  outputs: 


d 


o 


2  A 


sub 


(5.31) 
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We  emphasize  that  optimal  processing  within  clusters  is  usually 
realized  with  a  substantial  increase  in  instrumentation  complexity, 
especially  if  is  large.  One  ooviously  minimizes  this  complexity 
with  2  sensor  clusters  (i.e.,  =  2,  in  Fig.  4.2). 

Most  practical  instrumentations  will  depend  on  the  actual  spatial 
structure  of  the  noise  field,  which  usually  varies  with  towed  array 
conditions  (e.g.,  tow  vessel  speed).  Therefore,  adaptive 
instrumentations  are  anticipated  for  either  optimal  or  conventional 
processing  within  clusters,  in  conjunction  with  conventional 
beamforming  followed  by  the  best  postbeamformer  filter.  Modified 
Ecxart  filtering,  with  conventional  processing  within  clusters,  avoids 
excessive  instrumentation  complexity  associated  with  optimally 
processed  clusters  y.see  Fig.  4.2)  since  it  demands  only  a  single 
adaptive  frequency  scaling  at  the  conventional  beamformer  output  for 
homogeneous  noise.  From  Chapter  4  (see  Fig.  4.12),  the  general  post- 
beamforiner  filter  construction,  in  the  latter  case,  requires  M1  -  1 
bandpass  filter  parameters  ^center  frequencies)  with  sensors  per 
cluster.  For  sufficiently  large  M.,  we  shall  demonstrate  that  a 
single  ideal  bandpass  filter  with  simple  Ecxart  filtering,  at  a 
conventional  beam  output,  recovers  a  large  fraction  of  the  detection 
loss  incurred  when  only  simple  tcxart  filtering  is  used  at  a  beam 
output. 
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5.4  SuDarray  Size  and  Separation 

The  issue  of  optimum  subarray  (group)  size  and  subarray 
separation  with  a  limited  number  of  available  sensors  remains  to  be 
examined.  Fig.  5.4  shows  normalized  detection  indices  for  the 
detector  with  conventional  beamforming  on  optimally  or  conventionally 
processed  clustered  sensors  as  a  function  of  =  M/N^  sensors  per 
cluster,  using  overall  array  gain  expressions  from  (5.27)  and  (5.30), 
with  detection  indices  in  (5.28)  and  (5.31).  There  are  N  =  Tw/2tt  =  50 
positive  frequencies  in  the  signal  band  and  wa/Uc  =  2-n.  Signal-to- 
noise  ratio  is  again  assumed  uniform  in  the  signal  bandwidth  as 
specified  in  (5.4).  Hence,  signal -to-noise  ratio  and  the  M  =  N,M^ 
fixed  number  of  available  sensors  are  conveniently  normalized  out  of 
the  results.  In  all  curves  snown,  &  1  =  e  (M,  -  1)a  <  1.  As 

K  K  X 

one  snould  expect,  the  detection  index  with  optimal  processing  on 

9 

clustered  sensors  (dQ  )  is  always  larger  than  the  detection  index  with 

2 

conventional  processing  within  clusters  (dc  |modified  Eckart)-  We  see 

2  2 

tnat  all  curves  for  both  d  '  and  d  |  .  increase 

o  c  (modified  EcKart 

monotonical ly  witn  M.,  with  the  largest  improvement  variation 
occurring  from  one  to  two  sensors.  From  tne  monotonic  behavior,  one 
concludes  that  detection  performance  is  best  wnen  is  as  large  as 
is  practically  possible.  Given  a  fixed  M  number  of  available  sensors, 
M  =  ,  one  obviously  maximizes  by  minimizing  the  number  of 

clusters,  N, .  If  is  taken  large  enough  so  that  1  =  (M^  -  1)a 
exceeds  the  noise  correlation  length  (1  »  1/6^),  then  the  results 
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of  Chapter  4  indicate  that  conventional  processing  on  the  clustered 
sensors  is  near  optimal.  To  achieve  the  large  br1  condition 
obviously  requires  either  many  sensors  per  cluster  or  greater  sensor 
spacings  within  the  cluster;  with  the  parameters  of  Fig.  5.4,  the 
required  M^  is  greater  than  16. 

Note  that  performance  is  always  better  when  available  sensors  are 
clustered  instead  of  simply  spread  out  uniformly  over  the  observation 
interval  if  the  latter  arrangement  maxes  the  noise  uncorrelated  from 
sensor  to  sensor.  Relative  wideband  performance,  with  uncorrelated 
noise  (i.e.,  =  1  in  Fig.  5.4),  is  10  log  (N)  =  17  decibels  for 

N  =  50  positive  frequencies.  The  detection  index,  in  this  case,  is 
substantially  below  any  of  the  detection  indices  obtained  with  tightly 
grouped  sensors. 

Since  good  detection  calls  for  large  M^,  hence  small  N-^, 
whereas  localization  demands  >  3  [39]  -  [42],  a  logical  choice  is 
1^  _  3  With  M^  =  M/3  uniformly  spaced  sensors  per  conventionally 
processed  cluster.  Specifically,  one  would  locate  a  cluster  at  the 
observation  interval  end  points,  with  a  cluster  at  the  interval 
midpoint.  It  should  be  clear  that  regular  cluster  spacing  is  not 
necessary.  The  important  restriction  is  that  the  noise  is 
uncorrelated  between  cluster  outputs. 

The  only  argument  for  doing  something  other  than  simple 
conventional  processing  on  clustered  sensors  is  to  recover  the 
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difference  between  dQ2  and  dJmodified  £ckart  in  Fig.  5.4.  Again,  one 
accomplishes  this  performance  recovery  with  optimal  processing  within 
clusters.  If  the  optimal  cluster  approach  is  elected,  then  2  sensor 
clusters  offer  a  relatively  simple  optimal  instrumentation  over  the 
>>  1  optimal  instrumentation.  The  complexity  advantage  over 
selecting  >  2  is  obvious.  Physically,  it  is  probably  easier  to 
instrument  many  2  sensor  clusters  than  fewer  clusters  with  many 
sensors  per  cluster.  The  obvious  disadvantage  with  =  2  is  the 
tradeoff  of  approximately  6  decibels  of  performance  in  this  example 
relative  to  the  optimal  performance  with  =  16. 

To  clarify  the  functional  dependence  of  the  detection  indices  on 

Mp  Fig.  5.4d  presents  the  results  for  a  =  0.001  on  log-log  scales. 

We  see  that  conventional  beamforming  on  optimally  processed  clusters 

increases  by  approximately  3  decibels  witn  each  doubling  of  sensors 

per  cluster  for  2  <  M.  <  8.  For  larger  M^,  the  optimal  clustering 

curve  saturates  approximately  at  70  decibels,  whicn,  after  careful 

2  2  2  2 

examination  of  (5.28),  is  seen  to  be  approximately  N/a  n  (50/a  it  ). 

The  curve  for  conventional  detection  performance  with  the 
modified  Ecxart  filter  increases  with  a  slope  of  approximately 
one-half.  Performance,  in  this  case,  increases  as  while 
instrumentation  complexity  increases  proportionally  with  -  1.  In 
general,  we  require  -  1  bandpass  filter  center  frequencies  to 
construct  tne  postbeamformer  filter  (e.g.,  see  Fig.  4.12).  From  this, 
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we  infer  that,  given  fixed  M,  tne  performance /complexity  ratio  for 
conventionally  processed  clustered  sensors  increases  as  1/nTmi  when 


is  relatively  smal 1 . 


Since  all  subarrays  have  tne  same  array  gain,  one  only  needs  a 

single  filter  (introduced  after  tne  final  beamformer).  The  complexity 

one  is  then  concerned  with  is  that  of  frequency  dependence  of  the 

subarray  gain.  However,  the  problem  associated  with  complexity  of  the 

postbeamformer  filter  may  be  less  severe  than  it  seems.  Fig-  4.12 

suggests  that  array  gain  fluctuations  within  one  2n  interval  of 

w  a/U  become  smaller  as  the  number  of  sensors,  M, ,  increases. 

K  c 

If  this  is  true,  tnen  tne  postbeamformer  filter  only  needs  to 
eliminate  the  normalized  frequency  neighoorhoods  of  / Uc  =  2nn, 

n  =  0,  1. 

To  demonstrate  tnis,  we  nave  computed  the  detection  indices  in 

(5.28),  (5.31),  and  (2.77),  using  array  gains  from  (5.27)  and  (5.30), 

with  an  ideal  symmetrical  Dandpass  filter  of  unit  amplitude,  centered 

at  co  a/U  =  it,  applied  to  the  output  of  the  conventional 
k  c 

beamformer.  Detection  indices  are  plotted  in  Fig.  5.5  on  the  vertical 
axis  versus  percent  of  filter  bandwidth  to  total  processing  bandwidth 

wa/Uc  =  2tt. 

All  curves  are  normalized  by  the  square  of  the  number  of  total 
sensors  times  the  assumed  constant  signal-to-noise  ratio. 
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The  rummer,  M^,  of  sensors  used  is  25,  and  the  noise  fractional 
bandwidth,  in  each  case,  is  0.01. 

The  following  observations  are  plain.  The  detector  with 
conventional  beainforrning  on  optimally  processed  clustered  sensors 
achieves  the  highest  performance  value.  The  penalty  paid  for  using 
simple  Eckart  filtering  following  conventional  beamforming  is  a 
33-decibel  loss  in  detection  performance  relative  to  conventional 
detector  performance  with  tne  modified  Eckart  filter  when  the  percent 
bandwidth  of  the  bandpass  filter  is  sufficiently  large  (i.e.,  96 
percent)  to  include  aliased  frequencies  in  the  neighborhoods  of 
co,a/U  =  0  and  2n  (e.g.,  see  Fig.  4.12).  Aliased  frequencies  are 
clearly  eliminated  as  the  width  of  tne  bandpass  filter  relative  to  2w 
decreases.  Specifically,  the  uncomplicated  bandpass  filter  brings  the 
detector  with  conventional  beamforming  followed  by  simple  Eckart 
filtering  to  approximately  2.5  decibels  below  the  performance  with 
modified  Eckart  filtering  when  the  percent  of  total  bandwidth  is  50. 

At  the  percent  of  total  bandwidth  of  50,  optimal  detection  index  and 
detection  index  of  the  conventional  detector  with  modified  Eckart 
filtering  essentially  achieve  tneir  maximum  values  for  this  example. 

In  contrast,  performance  of  the  conventional  detector  with  simple 
Ecxart  filtering  begins  to  degrade  as  the  filter  bandwidth  widens  to 
admit  frequencies  in  aliased  neighborhoods ,  near  u>  a/U  =  0  and 

K  C 

wa/Uc  =  2ir .  Finally,  as  the  bandpass  filter  width  relative  to  2ir 
decreases  to  0.16,  conventional  detection  indices  converge  on  each 
other,  as  they  must  in  the  extreme  of  very  narrowband  processing. 
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Changing  conditions  naturally  lead  to  adaptive  instrumentations. 
One  thus  needs  only  a  single  adaptive  modified  Eckart  filter  with 
conventional  processing  on  clustered  sensors  witn  homogeneous  noise. 

On  the  other  hand,  optimal  processing  within  clusters  requires 
adaptive  processing  on  N-^M^  sensors.  Practically  speaKing,  a 
modified  Eckart  filtering  adaptive  instrumentation  is  probably  best 
implemented,  at  least  with  the  Corcos  flow  noise  model,  by 
constructing  bandpass  filters  at  a  beam  output  requiring  knowledge  of 
only  the  noise  speed,  U^.  For  a  large  number  of  sensors,  a  single 
bandpass  filter  centered  symmetrically  on  wkA/Uc  =  »  may  suffice. 

Again,  one  needs  only  noise  speed,  Uc,  to  construct  such  a  filter. 

There  is  one  otiier  degree  of  freedom  —  the  choice  of  the 
interval  a.  As  one  increases  a,  one  increases  wa/Uc.  The  evidence 
in  Chapter  4  indicates  that  increases  in  wa/Uc  Deyond  2u  have  little 

influence  on  dc2|  inodified  Eckart  (see  Fi9-  4-6b>  because  of  the  near 
periodicity  of  G.  ^  with  w.  a/U  .  By  increasing  a,  one  can  make 

K  C  K  C 

Sl  >  1,  thus  ensuring  that  Gkc  =  Gk3,  as  in  Fig.  4.4,  and 

2  2 

recapturing  all  of  tne  differences  between  dQ  and  dc  | mod i f i ed  Eckart 
in  Fig.  5.4.  This  would  lead  to  the  following  prescription  for  array 
design:  divide  the  sensors  into  three  clusters.  Choose  the  sensor 
separation  in  each  cluster  such  tnat  &kl  =  1  at  the  lowest  frequency 
and  highest  noise  speed.  Then  all  one  needs  is  a  simple  conventional 
beamformer  with  one  modified  Eckart  filter  at  the  output  to  achieve 
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near  optimal  performance.  This  near  optimal  performance  is  slightly 
poorer  tnan  the  performance  obtainable  with  tighter  spacing  of 
sensors,  but  the  implementation  is  much  simpler.  The  only  negative 
factor  is  the  complexity  of  the  modified  Eckart  filter,  which,  as  we 
have  demonstrated,  may  be  avoided  if  one  uses  a  simple  bandpass  filter 
at  the  beamforiner  output. 
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N,  EXTENDED  N,  CHANNELS 
SENSORS 

L  =  N,^+  (N.  -  I)  D 


Fiq.  5.1.  Conventional  Detector  Structure  With  Conventional 
Beamforming  on  N'  Space-Continuous  Extended  Sensors 
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Fig.  5.2.  Conventional  Detection  Enhancement  With  Conventional 
Beamforming  on  Extended  Sensors 
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N,  =  M  M,  (Clusters) 


L  =  Nt-£+  (Ni  I)  D 


Fig.  5.3.  Conventional  Detector  Structure  With  Conventional 
Beamforming  on  N,  Subarrays  of  Discrete  Sensors 
jf  =  (Mi  —  1)A,  M  =  NiMi  Sensors 
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a. 


c. 


b.  d 


Fig  5.4.  Conventional  Wideband  Detection  Indices 
for  Optimal  and  Conventional  Processing  Within  Clusters, 
w a  Ui  =  2tt,  or  —  fik  \v.  M  —  N i  M  i  and  Sk  K  =  S  I  (Co nstant) ( 
k  =  1.2,— N  Positive  Frequencies 
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Fig.  5.5.  Wideband  Detection  Indices  Versus  Percent  of  Total  Band, 
Aw  IT  =  2 77  With  Symmetric  Bandpass  Filter  Centered 
at  oaA  IT  =  7 r.  Constant  Signal  to  Noise  in 
[0,  w,]  and  N  =  50  Positive  Frequencies 
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CHAPTER  SIX 


CONCLUSIONS 


6.1  Summary  of  Results 

A  systematic  approach  has  been  developed  as  a  means  of 
characterizing  optimal  and  conventional  space-time  detection  of  a 
target  in  a  flow-noise-dominated  environment.  We  derived  the  optimum 
structure  for  a  stationary  Gaussian  signal  and  Gaussian  spatially 
"colored"  flow  noise  characterized  as  both  stationary  in  time  and 
homogeneous  in  space.  Signal  and  noise  were  assumed  independent  of 
each  other.  We  imposed  the  additional  constraint  that  the  signal 
propagated  as  a  plane  wave  from  a  Known  direction.  The  signal 
wavefront  was  assumed  to  travel  across  the  array  in  a  time  short 
compared  with  the  observation  interval  T,  with  signal  statistics 
satisfying  Tw  >>  1.  If,  in  addition,  the  noise  traveled  its 
correlation  length  in  a  time  short  compared  to  T  and  T  was  large 
compared  with  the  noise  correlation  length,  the  space-time  detection 
problem  became  a  problem  of  spatial  processing  at  independent  signal 
frequencies.  This  then  enabled  us  to  focus  our  attention  on  the 
spatial  aspects  of  the  optimal  processor  and  its  performance  at 
independent  signal  frequencies  under  tne  assumption  of  small  signal-to- 
noise  ratio  at  each  sensor. 
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The  wideband  processor  and  its  performance  followed  from 
appropriately  weighted  sums  of  corresponding  narrowband  components  at 
the  independent  signal  frequencies. 

We  formulated  exact  expressions  for  the  likelihood  ratio  test  for 
discrete  and  continuous  oDservations .  In  each  case,  tne  optimal 
instrumentation  and  its  performance  were  related  to  the  temporal  and 
spatial  properties  of  signal  and  Butterworth  rational  flow  noise 
model.  We  then  analyzed,  in  some  detail,  optimal  performance  and  also 
tne  best  possible  performance  advantage  attainable  witn  optimal 
processing  over  conventional  processing  (detector  with  conventional 
beamforming)  in  the  modeled  flow  noise  environment. 

In  passing  from  discrete  to  continuous  but  finite  length 
observations,  we  establisned  that  the  uncorrelated  noise  component 
formally  has  a  spectrum  level  of  zero.  Hence,  only  wavenumber-colored 
flow  noise  remained  to  limit  detection  when  one  had  continuous 
observations.  The  optimal  processor  structure  was  found  by  solving  a 
nontrivial  Fredholm  integral  equation  of  the  first  kind  encountered  in 
the  binary  detection  problem.  It  was  shown  to  consist  of  a 
conventional  beamformer  in  addition  to  appropriately  weignted 
end-point  sampling  sensors  located  at  each  end  of  the  finite 
observation  interval. 
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We  found  that  the  product  of  the  spatial  oandwidtn  and  array 
length  ( L )  is  a  very  important  spatial  parameter.  Array  gain  with 
conventional  beamforming  alone  was  shown  to  be  at  most  1.8  decibels 
below  optimal  with  the  first-order  Butterworth  model  if  the  array 
length  exceeds  two  noise  correlation  lengths.  This  very  important 
result  states  that  once  the  large  condition  holds,  nothing  more 
can  be  done  with  any  other  processor  to  significantly  improve  array 
gain  over  what  has  already  been  accomplisned  by  the  conventional 
beamformer  alone.  The  detection  problem  thus  must  be  reduced  to  that 
of  near  spatially  white  noise  in  a  band  of  order  1/L,  over  which  the 
response  of  the  conventional  beamformer  is  significant. 

We  indicated  tnat  substantial  array  gain  improvement  is  possible, 
over  conventional  beamforining  alone,  if  one  implements  end-point 
sampling  sensors  in  addition  to  conventional  oeamforming.  However, 
end-point  sensors  require  a  relatively  complex  instrumentation,  and 
they  significantly  improve  performance  only  when  the  array  length  is 
much  less  than  the  noise  correlation  length.  Sucii  a  condition  is 
proDably  best  avoided  in  many  towed  array  applications  by  simply 
extending  the  observation  interval  (array  length)  beyond  the  noise 
correlation  distance. 
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Even  the  idea  of  conventional  steering  ^beamforming)  is  likely  to 
be  impractical  over  any  space-continuous  length.  Therefore,  the 
formulation  of  space-continuous  arrays  with  any  steering  at  all 
probably  represents  an  unrealistic  idealization.  Its  purpose  is 
merely  to  establish  absolute  performance  bounds.  Such  bounds  are  used 
as  benchmarKS  for  comparisons  with  the  performance  of  other  array 
geometry.  Of  course  tnese  performance  bounds  can  never  be  reached 
with  any  other  array  geometry  whatsoever. 

A  key  development  step  noted  that  array  gain,  a  frequently  used 
measure  of  performance,  is  independent  of  a  linear  transfer  function 
connected  to  the  oeamformer  output,  but  wideband  detection  or  tracker 
performance  is  strongly  dependent  on  linear  filtering  applied  to  a 
conventional  beamformer  output.  For  flow  noise,  usually  concentrated 
at  high  wavenumbers,  optimal  and  conventional  detectors  can, 
therefore,  have  dramatically  different  wideband  detection  indices, 
even  when  conventional  array  gain  is  itself  nearly  optimal. 

We  nave  shown  that  much  of  the  performance  loss  due  to  the  use  of 
conventional  beamforming  can  be  recovered  by  optimal  postbeamformer 
scalar  frequency  filtering.  For  continuous  observation  intervals, 
long  compared  with  noise  correlation  length,  the  best  postbeamformer 
filter  can  reduce  the  difference  between  optimal  and  conventional 
( beamf ormer )  detector  performance  to  a  negligible  amount.  For  the 
Corcos  flow  noise  model,  the  best  postbeamformer  filter 
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instrumentation,  in  this  case,  is  a  simple  modification  of  the  basic 
EcKart  filter,  consisting  of  relatively  simple  pre-emphasis  of  high 
signal  frequencies.  A  remarkable  property  of  this  best  postbeamformer 
filter,  from  Appendix  A,  is  its  optimality  for  a  large  class  of  noise 
wavenumber  spectra  of  the  form 


where  v  and  e  are  linear  functions  of  frequency  &>  .  We  note 

K  K  K 

that  this  class  of  spectra  includes  the  simple  single  pole  Butterworth 
wavenumber  spectrum  at  one  extreme  ana  an  infinite  pole  Gaussian 
wavenumoer  spectrum  at  the  otner. 

In  Chapter  4,  we  examined  tne  likelihood  ratio  test  and  its 
performance  with  arrays  of  uniformly  spaced  sensors  in  the  absence  of 
uncorrelated  noise.  We  justified  omitting  uncorrelated  noise  with 
arrays  by  arguing  that  performance  predictions  without  the 
uncorrelated  noise  represent  the  upper  bound  performance  advantage 
achievable  with  optimal  processing  over  conventional  processing.  Even 
the  smallest  amount  of  uncorrelated  noise  tends  to  make  the  noise 
wavenumber  spectrum  flatter,  thus  bringing  the  conventional  beamformer 
closer  to  optimal . 


159 


TR  7997 


We  demonstrated  that  spatial  aliasing  is  ideally  all  that  happens 
when  one  goes  from  continuous  to  equally  spaced  discrete  observations. 
Aliasing  causes  high  wavenumDer  noise  components  to  fold  down  onto 
lower  acoustic  wavenumbers.  This  necessarily  increases  the  noise 
level  in  the  neighborhood  of  the  acoustic  signal  wavenumber  and  hence 
degrades  signal-to-noise  ratio  relative  to  continuous  observations. 

If  one  utilizes  the  common  practice  of  spacing  sensors  at  half 
acoustic  wavelengths,  which  is  often  much  larger  than  the  flow  noise 
correlation  length,  then  the  sampled  noise  wavenumber  spectrum  becomes 
essentially  flat.  Conventional  beamforming  is  nearly  optimal  at  this 
spacing,  but  beam  output  signal-to-noise  ratio  is  extremely  poor  since 
the  level  of  the  sampled  noise  spectrum  is  very  much  larger  than  the 
noise  level  with  continuous  ODservat ions . 

The  best  performance  attainable  with  continuous  observation  is 
nearly  achieved  witn  arrays  of  extremely  tight  sensor  spacings.  Such 
spacings  must  be  smaller  than  flow  noise  wavelengths,  which  typically 
have  a  length  scale  in  incnes.  One  avoids  spatial  aliasing  if,  and 
only  if,  spatial  observations  are  continuous.  Once  we  have  chosen  to 
discretize  the  observation  interval,  detection  performance  is  lost. 

No  amount  of  processing  applied  subsequently  can  recover  this  lost 
performance . 
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Aliasing  of  high  wavenumber  noise  is  related  directly  to 
variations  of  array  gain  with  frequency  when  the  noise  is  correlated 
between  sensors.  The  minima  of  these  variations  occur  at  frequencies 
where  sensor  separations  are  such  that  the  noise  components  add 
coherently.  The  maxima  occur  at  frequencies  where  noise  components  at 
adjacent  sensors  are  out  of  phase  and,  therefore,  tend  to  cancel 
naturally.  Noise  components  at  the  former  frequencies  can  dominate 
the  conventional  beamformer  output.  The  optimal  detector,  therefore, 
outperforms  the  detector  with  conventional  beamforming  and  simple 
Eckart  filtering  by  large  amounts.  If  the  undesirable  frequencies  are 
eliminated,  the  best  postbeamformer  filter,  optimized  for  its  actual 
input  signal-to-noise  ratio,  achieves  dramatically  improved  wideband 
performance  of  tne  detector  using  conventional  oeamforming.  This 
performance  becomes  nearly  optimal  when  conventional  array  gain  is 
nearly  optimal  at  all  signal  frequencies. 

This  treatment  of  arrays  suggests  that  the  use  of  a  relatively 
simple  modification  of  the  Eckart  filter  can  very  substantially 
enhance  the  performance  of  a  detector  using  conventional  beamforming, 
making  it  more  nearly  competitive  with  the  optimal  detector.  Major 
gains  can,  of  course,  be  made  only  when  a  reasonably  wide  band  of 
frequencies  is  processed  and  the  noise  has  sufficient  spatial 
structure  (i.e.,  small  fractional  bandwidth)  so  that  optimal  and 
conventional  array  gains  vary  significantly  with  frequency.  The  flow 
noise  associated  with  towed  arrays  has  tne  general  form  of  the  model 
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used  in  tnis  study.  While  its  fractional  bandwidth  need  not  be  as  low 
as  the  values  considered,  it  would  probably  be  sufficiently 
concentrated  in  wavenumber  space  to  make  the  modified  Eckai t  filter 
attractive  for  wideband  detection. 

A  clear  drawback  of  the  modified  Eckart  filter  with  arrays  is,  of 
course,  its  dependence  on  the  (conventional)  array  gain  and  hence  on 
the  spatial  structure  of  tne  noise  field.  This  structure  is  not 
likely  to  be  known  in  advance  so  that  some  form  of  adaptation  is 
required.  However,  one  needs  only  a  single  adaptive  filter  at  a  beam 
output.  For  arrays,  this  suggests  a  significant  reduction  in 
processing  requirements  compared  with  the  set  of  M  filters  usually 
required  for  optimal  processing  of  data  from  M  sensors.  For  flow 
noise  encountered  by  towed  arrays,  the  primary  unknown  parameter  would 
probaDly  be  convective  noise  speed,  U^,  a  quantity  closely  related 
to  tow  speed.  Adaptation  could,  therefore,  be  a  relatively  simple 
matter. 

In  Chapter  5,  we  dealt  with  criteria  that  affect  the  choice  of 
geometry  in  practice.  With  continuous  observations,  the  primary 
constraint  is  that  the  introduction  of  continuous  phase  shifts  is 
impractical  so  that  one  must  work  with  sensor  lengths  that  are  small 
compared  with  acoustic  wavelengths. 


16? 


TR  7997 


The  impractical i ty  of  continuous  phase  shifts  led  us  to  consider 
a  detector  with  conventional  beamforming  on  space-continuous  subarrays 
(extended  sensors)  wnere  steering  within  a  sensor  is  abandoned. 

Our  results  with  the  simple  Butterworth  model  suggest  that 
conventional  Deamforming  on  a  fully  populated  (filled)  array  of 
individual  equally  spaced  extended  sensors,  with  lengths  sufficiently 
small  compared  with  signal  wavelengtns,  nearly  achieves  the  best 
possible  array  gain  performance  if  the  array  length  is  sufficiently 
larger  than  tne  noise  correlation  length.  One  then  nearly  achieves 
best  wideband  detection  index  if,  in  addition,  the  array  length 
exceeds  trie  noise  correlation  length  at  tne  lowest  processed  signal 
frequency. 

An  alternate  array  geometry  uses,  in  general,  irregularly  spaced 
extended  sensors  so  that  the  noise  is  uncorrelated  between  sensor 
outputs.  Then  conventional  beamforming  on  extended  sensor  outputs 
followed  Dy  modified  Eckart  filtering  can  come  close  to  optimal 
wideband  performance,  provided  that  each  extended  sensor  length  is 
less  than  one-fourth  of  the  snallest  acoustic  signal  wavelengtn  along 
the  array.  Processed  signal  oandwidtn  is  thus  kept  wnol ly  inside  the 
frequency  interval  where  the  extended  sensor  response  to  the  signal  is 
near  unity.  Extended  sensors  are  effectively  steered  in  this  case. 

If  eacn  extended  sensor  length  also  exceeds  the  noise  correlation 
length  at  the  lowest  signal  frequency  and,  in  addition,  the  unused 
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aperture  between  extended  sensors  is  a  small  fraction  of  the  entire 
observation  interval,  then  conventional  beamforming  on  equally  spaced 
extended  sensors,  in  conjunction  with  best  postbeamf ormer  filtering? 
nearly  achieves  the  best  detection  index  attainable  witn  continuous 
observations  on  the  entire  observation  interval. 

For  this  special  case  of  extended  sensor  lengths,  long  compared 
with  the  (Butterworth)  noise  correlation  length  and  simultaneously 
short  compared  with  signal  wavelengths,  the  modification  to  the  simple 
Ecxart  postbeamf  oriner  filter,  analogous  to  full  array  processing,  is  a 
simple  voltage  gain  increasing  linearly  with  frequency.  Again,  this 
postbeamf  oriner  filter  remains  optimal  as  the  input  noise  wavenumber 
spectra  ranges  over  tne  large  class  of  input  noise  spectra  mentioned 
earlier  in  connection  witn  full  array  processing. 

We  considered  the  important  sonar  issue  of  the  constraint  of  a 
limited  number  of  point  sensors  in  a  given  observation  interval.  The 
number  must  often  be  small  enough  so  that  uniform  spacing  over  the 
available  aperture  results  in  flow  noise  being  independent  from  the 
sensor.  Otherwise,  one  can  nearly  achieve  the  optimal  performance 
with  conventional  beamforming  on  a  filled  array  of  tightly  but  equally 
spaced  sensors  if  the  array  length  exceeds  the  noise  correlation 
length.  By  using  groups  ^clusters)  of  tightly  spaced  sensors,  one  can 
get  a  much  higher  array  gain  at  certain  frequencies  and  hence  achieve 
greatly  improved  detector  performance  Meitner  optimal  or  conventional ) . 
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Again,  the  best  postbeamforiner  frequency  filter  must  be  applied  to  the 
conventional  beamformer  output. 

With  conventional  beamforming  on  subarrays  (clusters)  of  tightly 
grouped  discrete  sensors,  conventional  steering  and  optimal  processing 
within  clusters  are  potentially  possiDle.  Thus,  there  is  no 
requirement  that  the  subarray  lengtn  exceed  the  noise  correlation 
length  in  order  to  achieve  optimal  array  gain  at  tne  cluster  level. 

The  subarrays  are  separated  so  that  the  noise  is  uncorrelated 
between  them.  Again,  the  best  postbeamforiner  frequency  filter,  which 
incorporates  the  common  subaperture  array  gain  and  simple  Eckart 
filtering,  recovers  much  of  the  wideoand  detection  loss  incurred  by 
simple  Eckart  filtering  following  conventional  beamforming  on 
subarrays.  Substantial  recovery  of  detection  performance  is  possible 
at  the  beam  output  when  a  reasonably  wide  oand  of  frequencies  is 
processed  so  that  there  is  a  large  variation  in  subarray  gain  over 
frequency. 

We  noted  that  detection  performance  with  a  fixed  numoer  of 
sensors  is  maximized  when  the  number  of  sensors  within  a  cluster  is 
maximum.  Good  detection  calls  for  a  large  number  of  sensors  in  a 
cluster,  whereas  localization  criteria  demand  a  minimum  of  three 
clusters.  The  logical  choice  is  one  cluster  at  each  end  of  the 
observation  interval  with  the  remaining  cluster  at  the  midpoint  of  the 
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observation  interval.  Cluster  separation,  however,  need  not  be 
uniform.  The  important  restriction  is  that  the  noise  is  uncorrelated 
between  cluster  outputs.  The  suggested  array  construction  of 
clustered  sensors  with  conventional  steering  within  clusters  thus 
achieves  greatly  improved  detection  and  localization  performance  over 
conventional  beamforming  on  available  sensors  spread  over  the 
observation  interval.  This  latter  case,  with  flow  noise  independent 
from  sensor  to  sensor,  yields  an  array  gain  of  only  >1,  the  total 
number  of  sensors  used. 

Optimal  processing  within  a  cluster  always  improves  detection 
performance  over  simpler  conventional  steering  within  clusters, 
especially  if  the  cluster  length  is  much  less  than  the  noise 
correlation  length  at  most  signal  frequencies.  However,  there  is  a 
substantial  cost  in  complexity,  especially  if  tne  number  of  sensors 
per  cluster  is  greater  than  2.  Thus,  an  efficient  detector  structure 
uses  conventional  beamforming  on  outputs  of  optimally  processed  2 
sensor  clusters  distributed  in  the  observation  interval  in  any  manner 
whatever,  with  noise  uncorrelated  between  cluster  outputs.  The  best 
postbeamf ormer  filter,  in  this  case,  consists  solely  of  the  simple 
Ecxart  filter  chosen  from  signal  and  noise  spectra  at  each  sensor. 

One  can  exercise  the  sensor  separation  degree  of  freedom.  If  one 
chooses  uniform  sensor  separation  wic'nin  each  cluster,  such  that  the 
cluster  length  exceeds  the  noise  correlation  length  at  the  lowest 
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signal  frequency  and  fastest  tow  speed,  then  all  one  needs  is  a  simple 
conventional  beamformer  on  outputs,  witn  one  modified  tcKart  filter  at 
the  beamformer  output,  to  achieve  near  optimal  wideband  detection 
i ndex . 

6.2  Concluding  Remarks 

We  modeled  flow  noise  as  a  single  pole  rational  Butterworth 
spectrum.  It,  instead,  we  used  a  rational  spectrum  with  higher  order 
poles,  then  the  solution  of  the  Fredholm  integral  equations 
encountered  in  our  detection  proolem  with  continuous  observations 
contains  Dirac  delta  functions  ana  their  derivatives  at  the  interval 
end  points.  It  follows  that  the  likelihood  ratio  test  then  requires 
data  samples  and  derivatives  of  the  data  at  the  interval  end  points. 
This  naturally  increases  the  complexity  and  sensitivity  of  the  test 
over  the  test  specified  with  tne  first-order  Butterworth  model. 
However,  asymptotic  optimality  of  the  conventional  array  gain,  with 
the  plane  wave  signal,  still  holds  approximately  if  the  array  length 
is  sufficiently  larger  than  the  noise  correlation  length.  The  plane 
wave  signal  model  of  Known  location  can  be  shown  to  satisfy 
approximately  the  frednolm  integral  equation  for  more  general  noise 
wavenumber  spectra,  provided  the  array  length  is  sufficiently  larger 
than  the  noise  correlation  length  [2,  19, 
pp.  204-207 j . 
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With  continuous  observations  and  rational  noise  spectral  models, 
one  easily  argues  that  the  likelihood  ratio  test  designed  with  assumed 
noise  spectral  forms  can  be  very  sensitive  to  the  functional  form  of 
actual  noise  spectra.  Finite  power  flow  noise,  considered  to  be  in 
the  class  of  rational  spectra,  can  oe  expressed  as  a  ratio  of 
polynomials,  with  difference  p  -  q  >  1  between  the  order  of  the 
denominator  and  numerator  polynomials.  If  one  restricts  processing  to 
the  conventional  beamformer  alone,  then  its  performance  will  probably 
not  vary  much  over  the  class  of  rational  noise  spectra  once  the  total 
array  length  sufficiently  exceeds  the  noise  correlation  length. 

However,  this  is  not  true  for  strictly  optimal  processors 
designed  on  the  basis  of  an  assumeJ  functional  form  for  the  flow 
noise.  Performance  with  the  optimum  processor,  requiring  derivatives 
of  the  end-point  data,  can  vary  drastically  over  the  class  of  rational 
spectra  if  the  actual  noise  spectra  deviate  from  the  assumed  form, 
particularly  if  the  order  of  the  assumed  functional  form  in  the 
optimal  design  exceeds  the  order  of  the  actual  noise  spectra.  This 
suggests  a  way  to  make  the  theoretical  optimal  processor  "robust 
(less  sensitive).  One  makes  the  optimal  processor  more  robust  by 
reducing  or  eliminating  the  order  of  tne  data  differentiation  at  the 
interval  end  points.  With  rational  noise  spectra,  one  achieves  this 
with  lowest  order  optimal  processor  designs,  such  as  the  first-order 
Butterworth  wavenumber  spectrum  used  in  this  study.  This  then 
suggests  that  in  any  practical  problem  where  actual  noise  spectra  are 
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seldom  known,  it  is  prooably  wise  to  instrument  optimal  processors 
assuming  tne  lowest  possible  (nontrivial)  order  functional  form  for 
the  noise  spectra. 

6.3  Suggestions  for  Future  Work 

In  almost  all  cases,  we  have  obtained  exact  solutions  to  our 
formulation  of  the  binary  detection  problem  in  flow  noise.  As  usual, 
there  remain  unresolved  issues,  a  partial  list  of  which  follows: 

1.  We  have  not  established  mathematically  (for  wavenumber 
spectra  more  general  than  the  first-order  Butterworth)  that 
conventional  beamforming  on  a  plane  wave  signal  with  arrays  of 
discrete  sensors  is  asymptotically  optimal  when  the  array  length  is 
sufficiently  large  compared  with  noise  correlation  length.  Our 
mathematical  analysis  and  numerical  computations  indicate  that 
conventional  beamforming  is  nearly  optimal  when  the  array  length 
exceeds  two  noise  ( Butterworth)  correlation  lengths.  If  one  can 
establish  this  for  general  noise  spectra,  one  can  assert  that  beams 
not  carrying  the  signal  cannot  be  used  to  significantly  reduce  noise 
on  target-carryi ng  beams  if  the  array  length  is  long  compared  to  the 
noise  correlation  length.  Under  these  circumstances,  only  the 
signal-carrying  beams  are  relevant  for  detection. 
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2.  We  have  not  considered  the  complete  sensitivity  problem, 
specifically  amplitude  and  phase  error  tolerance,  which  will  surely 
impact  any  of  our  predicted  performance  values.  This  is  a  basic  issue 
for  any  paper  study.  We  feel,  however,  that  consideration  of  this 
issue  will  serve  to  reinforce  the  argument  for  conventional 
beamforming  followed  by  modified  tcnart  postbeamforiner  filtering.  It 
frequently  happens  in  practice  that  such  suboptimum  instrumentations 
are  least  sensitive  to  error  tolerances.  In  contrast,  optimal 
instrumentations  are  usually  more  complex  to  build  than  the  simpler 
conventional  processor  and,  typically,  are  more  sensitive  to  error 
tolerances. 


170 


TR  7997 


APPENDIX  A 
DETECTION  INDEX  WITH 
CONTINUOUS  OBSERVATIONS 

The  space-continuous  form  of  the  conventional  detection  index  for 
detector  structure  in  Fig.  2.5,  with  frequency  scaling,  is  ODtained 
directly  from  (2.78) : 


K  =  1 


?  *  * 

S,  H  V,  VL  ht 

k  — <  -^k  — k  -k 


E 

k=l 


cj4  n* 

k  k 


itc'  ^  ilk  2 


( A.  1 ) 


Inserting  ,  the  conventional  detection  -index  with 

frequency  scaling  is 


d2  = 


E  K 

k=l  1 

2 

N 

E  K  4  \2 

k=l  1 

M  M 

E  E  *i’  Qii  vi  vi 

i=l  1=1 

2 

(A. 2) 


A  =  AiA1  , 
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where  i  and  1  denote  identical  uniform  sampling  of  interval  A  on  x  and  y, 

respectively,  in  [-L/2,  L/2].  Witn  N^Q.^  =  qk(x.,y-j),  a 

general  sampled  noise  covariance  function,  we  write  (A. 2)  as  follows: 


d2  = 


2  2 
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,  12  c 
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*1 
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(A. 3) 


Taking  the  limit  of  (A. 3)  as  A-,  A|  >  o,  with  continuity 
conditions  assumed  satisfied,  we  obtain 


(A. 4) 


where 
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For  sufficiently  large  T,  the  sum  over  k  can  be  approximated  by  an 
integral  yielding 


1  fw  1  1 

2 

L  I  c(w) 

1  Jo  \ 

2  S(io)dijj 
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|CU)|4 

Jq  1 

L 12 

f  C  +jvsy  -jvsx 

JJ  e  e  q(aj;x,y)dxdy 

-L 12 

2 

da) 

(A. 5) 


where  q^(x,y)  >  q  (u>;x,y),  >  C(w),  and  it  is  understood  that 

\>s  is  a  function  of  radian  frequency  «». 

Array  gain  for  tne  conventional  detector  structure  on 
space-continuous  [-L/2,  1/2]  is  obtained  easily  by  considering  only 
tne  ktn  component  of  (A. 4): 


A 


qk(x,y)dxdy 


(A. 6) 


Note  that  (A. 6)  is  not  affected  by  frequency  scaling,  C^.  Hence, 
narrowband  conventional  array  gain  is  not  influenced  by  frequency 
scaling  following  conventional  beamforming. 
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Let 


x  -  y  =  a 

(A. 7) 


x  +  y  =  y  ; 


then  the  unique  solution  for  x  and  y  is 


x 


x 


(A. 8) 


The  Jacobian  of  the  transformation  is  1/2.  It  follows  that 


-jv  J 

do  e  5  qK(j)(L+o) 
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(A. 9) 
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or,  the  equivalent. 


-jvsx 

e  qK(x-y)dxdy 


da 


(A. 10) 


Inserting  (A. 10)  into  (A. 6),  we  obtain 
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Letting  e.L  >>  1, 
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where  we  have  used  the  Wiener-Knintcmne  relationship  [37]  between  tne 
noise  spatial  covariance  function  qk(o)  and  its  wavenumber  spectrum 
Qn(vs).  Inserting  (A. 12)  into  (A. 11),  we  have  the  large 
conventional  array  gain  for  the  general  noise  wavenumber  spectrum 
Qn(v)  with  power  spectrum  1^: 


L  I 


’kc  -  q^T  • 


(A. 13) 


8kL  »  1 


For  the  first-order  Butterwortn  Homogeneous  covariance  function  in 
(2.34),  we  have,  from  (A. 10),  that 
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These  are  tabulated  integrals  |_38,  Gradshteyn  and  Ryzhik,  p.  196 
no.  3,  p.  198,  no.  6],  which  combine  to  yield 
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Inserting  (A. 15)  into  (A. 6),  we  obtain 
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or,  equivalently, 


(A. 17) 
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APPENDIX  6 
ARRAY  GAIN 

WITH  DISCRETE  OBSERVATIONS 


For  a  linear  array  of  M  regularly  spaced  sensors  i 
separated  by  a  distance. 


Tne  elements  of  the  first-order  Butterwortn  noise 
matrix  Qk  in  (4.6)  from  l 2. 34)  are 


-3k|(r-n)A  |  +  jvk(r-n)a 
6  > 

r,n  =  1,  2,  ....  M 


Optimal  Array  Gain 

Optimal  array  gain  from  (2.49)  is 


G. 


KO 


-  V  *  0  ^  V 

-  \  \  \ 


wnere  the  "steering  vector"  V^  is  defined  by  (2.22). 
and  (8.3),  we  obtain 
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From  (2.22) 
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Diagonal  elements  of  banded,  Herinitian  Q,  ^  contribute  terms 
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Conventional  Array  Gain 


Conventional  array  gain  is  ootained  from  the  general  expression 
in  (2.69): 
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M" 


KC 
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(8.6) 


Expanding  the  quadratic  form  in  the  denominator  of  (6.6),  using  Qk 
from  (4.6),  we  nave 
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An  equivalent  form  for  (6.7)  is 
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Inserting  (6.8)  into  (B.6),  we  arrive  at 
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A  more  useful  expression  for  (B.9)  takes  tne  form 
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Sums  (geometric  progressions)  in  the  denominator  of  (B.lOa)  have 
closed  forms  from  [38,  p.  1,  Equation  0.113],  whicn,  when  inserted 
into  (B.lOa),  yield 
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A  straightforward  calculation  yields  an  equivalent  expression, 


(B. 11a) 


Array  Gain  for  Conventional 
Beamforming  on  Subarrays 

In  physical  terms,  the  overall  array  gain  for  tne  clustered 
structure  is  the  signal-to-noise  ratio  at  tne  beamformer  output 
divided  by  the  signal-to-noise  ratio  at  each  sensor. 

Tne  array  gain  for  a  cluster  is  obtained  from  the  general 
expression  for  array  gain  in  (2.63): 
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The  matrix  filter  for  each  cluster  steers  to  wavenumber  vQ: 
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where 


VD  = 


cos  © 


D 


(3.14) 


Note,  \>q  does  not  necessarily  coincide  with  the  signal  wavenumber 
v$.  One  determines  array  gain  for  a  cluster  using 
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where  now  ^  •  x_/c  is  the  delay  at  point  x  in  the  cluster  relative  to 
the  origin  or  reference  point  within  a  cluster.  For  uniformly  spaced 
sensors,  assuming  that  the  first  sensor  within  a  cluster  is  the 
reference  point, 
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It  follows  from  (3.16)  and  the  properties  of  a  geometric  series  that 


(8.17) 


If  we  now  define 
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Note  this  is  not  the  usual  G^c  since  \>s  ^  in  (3.13). 

The  overall  array  gain  at  the  beamformer  output  is  obtained  in  a 
similar  manner  from  (2.63)  by  treating  the  clusters  as  individual 
sensors  with  noise  uncorrel ated  between  clusters.  With  beamformer 
matrix  filter 
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we  obtain 
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( B . 21 ) 


Since  the  noise  is  uncorrelatea  between  clusters,  we  have,  in  this 
case, 


Qk  =  I  (Identity  Matrix). 


(B.22) 
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1 


Analytically,  one  generates  overall  array  gain,  witn  correct  steering 
within  subarrays,  oy  setting  =  \>s-  Hence  ( 8 .23)  becomes 
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kD 


(B. 24) 


One  extends  (B.23)  to  array  gain  for  the  space-continuous 
suoaperture  in  fig.  5.1  in  the  manner  of  Appendix  A.  Taking  the  limit 
in  (B.23)  as  A  >  0,  with  continuity  conditions  assumed  satisfied  on 
fixed  1  and  setting  =  o  in  (B.17),  we  obtain 
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i 


represents  the  response  of  a  space-continuous  subaperture  (extended 
sensor)  with  generally  no  steering  witnin  tne  sensor  to  the  plane  wave 
signal  at  v$  unless  v$  =  0.  in  (B.13)  becomes  a  vector  of 
ones,  in  this  case,  and 
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kC|subc  ?Bk  .  ^iTV*  2e"Bk' 

\V  '(\S2)2  ><vkV)2 

(8.27) 

or,  equivalently. 


Equation  (B. 28)  is  analogous  to  the  array  gain  expression  in  (3.14)  for 
the  entire  array.  Tne  key  difference  in  ^3.28)  is  that  the  extended 
sensor  is  always  steered  to  broadside. 
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